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Abstract 

This paper describes a theory of transverse coherence properties of Undulator Radi- 
ation. Our study is of very practical relevance, because it yields specific predictions 
of Undulator Radiation cross-spectral density in various parts of the beamline. On 
the contrary, usual estimations of coherence properties assume that the undulator 
source is quasi-homogeneous, like thermal sources, and rely on the application of 
van Cittert-Zernike (VCZ) theorem, in its original or generalized form, for calculat- 
ing transverse coherence length in the far-field approximation. The VCZ theorem is 
derived in the frame of Statistical Optics using a number of restrictive assumptions: 
in particular, the quasi-homogeneous assumption is demonstrated to be inaccurate 
in many practical situations regarding undulator sources. We propose a technique 
to calculate the cross-spectral density from undulator sources in the most general 
case. Also, we find the region of applicability of the quasi-homogeneous model and 
we present an analytical expression for the cross-spectral density which is valid up 
to the exit of the undulator. For the case of more general undulator sources, simple 
formulas for the transverse coherence length, interpolated from numerical calcula- 
tions and suitable for beamline design applications are found. Finally, using a simple 
vertical slit, we show how transverse coherence properties of an X-ray beam can be 
manipulated to obtain a larger coherent spot-size on a sample. This invention was 
devised almost entirely on the basis of theoretical ideas developed throughout this 
paper. 
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1 Introduction 



In recent years, continuous evolution of third generation light sources has 
allowed dramatic increase of brilliance with respect to older designs, which 
has triggered a number of new techniques and experiments unthinkable before. 
Among the most exciting properties of today third generation facilities is the 
high flux of coherent X-rays provided. The availability of intense coherent X- 
ray beams has fostered the development of new coherence-based techniques like 
fluctuation correlation dynamics, phase imaging, coherent X-ray diffraction 
(CXD) and X-ray holography. In this context, understanding the evolution 
of transverse coherence properties of Synchrotron Radiation (SR) along the 
beam line is of fundamental importance. 

In general, when dealing with this problem, one should account for the fact 
that Synchrotron Radiation is a random statistical process. Therefore, the 
evolution of transverse coherence properties should be treated in terms of 
probabilistic statements: the shot noise in the electron beam causes fluctua- 
tions of the beam density which are random in time and space. As a result, 
the radiation produced by such a beam has random amplitudes and phases. 

Statistical Optics [1, 0, affords convenient tools to deal with fluctuating 
electromagnetic fields in an appropriate way. Among the most important quan- 
tities needed to describe coherent phenomena in the framework of Statistical 
Optics is the correlation function of the electric field. In any interference ex- 
periment one needs to know the system (second order) correlation function 
of the signal at a certain time and position with the signal at another time 
and position. Alternatively, and equivalently, one can describe the same ex- 
periment in frequency domain. In this case one is interested in the correlation 
function of the Fourier transform of the time domain signal at a certain fre- 
quency and position with the Fourier transform of the time domain signal 
at another frequency and position. The signal one is interested to study is, 
indeed, the Fourier transform of the original signal in time domain. In SR 
experiments the analysis in frequency domain is much more natural than that 
in the time domain. In fact, up-to-date detectors are limited to about 100 ps 
time resolution and they are by no means able to resolve a single X-ray pulse 
in time domain. They work, instead, by counting the number of photons at 
a certain frequency over an integration time longer than the radiation pulse. 
Therefore, in this paper we will deal with signals in the frequency domain and 
we will often refer to the "Fourier transform of the electric field" simply as 
"the field". 

For some particular experiment one may be interested in higher order correla- 
tion functions (for instance, in the correlation between the intensities) which, 
in general, must be calculated separately. In the particular case when the field 
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fluctuations can be described as a Gaussian process, the field is often said to 
obey Gaussian statistics. In this case, with the help of the Moment Theorem 
one can recover correlation functions of any order from the knowledge of 
the second order one: this constitutes a great simplification to the task of de- 
scribing coherence properties of light. A practical example of a field obeying 
Gaussian statistics is constituted by the case of polarized thermal light. This 
is more than a simple example: in fact, Statistical Optics has largely developed 
in connection with problems involving optical sources emitting thermal light 
like the sun, other stars, or incandescent lamps. As a consequence, Gaussian 
statistics is often taken for granted. Anyway, it is not a priori clear wether 
Synchrotron Radiation fields obey it or not; our analysis will show that Syn- 
chrotron Radiation is indeed a Gaussian random process. Therefore, as is also 
the case for polarized thermal light and any other signal obeying Gaussian 
statistics, when we deal with Synchrotron Radiation the basic quantity to 
consider is the second order correlation function of the field. Moreover, as al- 
ready discussed, in Synchrotron Radiation experiments it is natural to work in 
the space-frequency domain, so that we will focus, in particular, on the second 
order correlation function in the space-frequency domain. 

Besides obeying Gaussian statistics, polarized thermal light has two other spe- 
cific properties allowing simplifications of the theory: the first is stationarity 1 
and the second is quasi- homogeneity. Exactly as the property of Gaussian 
statistics, also stationarity and quasi-homogeneity of the source are usually 
taken for granted in Statistical Optics problems but, unlike it, they do not 
belong, in general, to Synchrotron Radiation fields. In fact, as we will show, 
Synchrotron radiation fields are intrinsically non-stationary and not always 
quasi-homogeneous. Nevertheless, up to now it has been a widespread prac- 
tice to assume that undulator sources are completely incoherent (i.e. homo- 
geneous) and to apply the well known van Cittert-Zernike (VCZ) theorem for 
calculating the degree of transverse coherence in the far-field approximation 
0]. Using the VCZ theorem, the electric field cross- correlation function in the 
far field is usually calculated (aside for a geometrical phase factor) as a Fourier 
transformation of the intensity distribution of the source, customarily located 
at the exit of the undulator. 

Although the VCZ theorem only deals with completely incoherent sources, 
there exists an analogous generalized version of it which allows to extend 
the treat the case of quasi-homogeneous sources as well. Actually there is 
no unambiguous choice of terminology in literature regarding the scope of 
the VCZ theorem. For instance, a very well-known textbook p reports of a 
" Zernike-propagation" equation dealing with any distance from the source. 
Also, sometimes the generalized VCZ theorem is referred to as Schell's 



Here we do not distinguish between different kind of stationarity because, under 
the assumption of a Gaussian process, these concepts simply coincide. 
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theorem (and also used in some paper [2j]). In this paper we will refer to the 
VCZ theorem and its generalized version only in the limit for a large distance 
from the source and for, respectively, homogeneous and quasi-homogeneous 
sources. However, irrespectively of different denominations, the fundamental 
fact holds, that once a cross-correlation function is known on a given source 
plane it can be propagated through the beamline at any distance from the 
source. It should be noted that, from this viewpoint, a source simply denotes 
an initial plane down the beamline from which the cross-correlation function 
is propagated further. Then, the position of the source down the beamline 
is suggested only on the ground of opportunity. On the contrary, when deal- 
ing with the VCZ theorem, the source must be (quasi)-homogeneous which 
explains the customary location at the exit of the undulator. 

In some cases, the VCZ theorem or its generalized version may provide a con- 
venient method for calculating the degree of transverse coherence in various 
parts of the beamline once the transverse coherence properties of the photon 
beam are specified at the exit of the undulator, that is at the source plane. In 
most SR applications though, such treatment is questionable. First, the source 
(even at the exit of the undulator) may not be quasi-homogeneous. Second, 
even for specific sets of problem parameters where the quasi-homogeneous 
model is accurate, the specification of the far-field zone depends not only on 
the electron beam sizes, but also on the electron beam divergencies (in both 
direction) and on the intrinsic divergence of the radiation connected with the 
undulator device. At the time being, widespread and a-critical use of the VCZ 
theorem and its generalization shows that there is no understanding of trans- 
verse coherence properties of X-ray beams in third generation Synchrotron 
Radiation sources. 

If, on the one hand, the definition of the far-zone and the possible non quasi- 
homogeneity of SR sources constitute serious problems in the description of 
the coherence properties of Synchrotron light, on the other hand the intrinsic 
non-stationarity of the SR process does not play a very important role. In par- 
ticular, as we will show, assumption of a minimal undulator bandwidth much 
larger than the characteristic inverse bunch duration (which is always verified 
in practice) allows to separate the correlation function in space-frequency do- 
main in the product of two functions. The first function is a spectral correlation 
describing correlation in frequency. The second function describes correlation 
in space and is well-known also in the case of stationary processes as the cross- 
spectral density of the process. Then, the cross-spectral density can be studied 
independently at any given frequency giving information on the spatial cor- 
relation of the field. Subsequently, the knowledge of the spectral correlation 
function brings back the full expression for the space-frequency correlation. 

In this paper we aim at the development of a theory of transverse coherence 
capable of providing very specific predictions, relevant to practice, regarding 
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the cross-spectral density of undulator radiation at various positions along 
the beam-line. A fully general study of undulator sources is not a trivial one. 
Difficulties arise when one tries to include simultaneously the effect of intrinsic 
divergence of the radiation due to the presence of the undulator, of electron 
beam size and electron beam divergence into the insertion device. The full 
problem, including all effects, poses an unsolvable analytical challenge, and 
numerical calculations are to be preferred. Generally, the cross-spectral density 
of the undulator radiation is controlled by nine physical parameters which 
model both the electron beam and the undulator: the horizontal and vertical 
geometrical emittances of the electron beam e X;y , the horizontal and vertical 
minimal betatron functions (3° y , the observation distance down the beamline 
z , the observation frequency u, the undulator resonant frequency u , the 
undulator length L w and the length of the undulator period X w . 

We will make a consistent use of dimensional analysis. Dimensional analysis 
of any problem, performed prior to analytical or numerical investigations, not 
only reduces the number of independent terms, but also allows one to clas- 
sify the grouping of dimensional variables in a way that is most suitable for 
subsequent study. The algorithm for calculating the cross-spectral density can 
be formulated as a relation between dimensionless quantities. After appro- 
priate normalization, the radiation cross-spectral density from an undulator 
device is described by six dimensionless quantities: the normalized emittances 
e Xj y = u e Xt y/c, the normalized betatron functions $ x>y = (3° /L w , the normal- 
ized observation distance z Q = z j L w and the normalized detuning parameter 
C = 27iN w (u—u )/lj , where N w = L w /\ w is the number of undulator periods. 

At some point in this work we will find it convenient to pose C = 0. In other 
words we will assume that parameters are tuned at perfect resonance. It is 
relevant to note that even under this simplifying assumptions, conditions for 
the undulator source to be quasi-homogeneous still include four parameters e XiV 
and P XjV . For storage rings that are in operation or planned in the Angstrom 
wavelength range, the parameter variation of e x ~ 10 — 10 3 , e y ~ 10 _1 — 10, 
Px,y ~ 10 _1 — 10 are possible: these include many practical situations in which 
the assumption of quasi-homogeneous sources and, therefore the (generalized) 
VCZ theorem, is not accurate. 

In this paper we will first deal with the most general case of non-homogeneous 
sources. In fact, from a practical viewpoint, it is important to determine the 
cross-spectral density as a function of e X;y , (3 X;y , C and z . Once a general 
expression for the cross-spectral density is found, it can be used as a basis 
for numerical calculations. A second goal of this work is to find the region of 
applicability of the quasi- homogeneous source model (i.e. of the generalized 
VCZ theorem) which will arise automatically from the dimensional analysis of 
the problem. Finally, we will derive analytical expressions for the cross-spectral 
density at C = in various parts of the beamline. 
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Results may also be obtained using numerical techniques alone, starting from 
the Lienard-Wiechert expressions for the electromagnetic field and applying 
the definition of the field correlation function without any analytical manipu- 
lation. Yet, computer codes can calculate properties for a given set of param- 
eters, but can hardly improve physical understanding, which is particularly 
important in the stage of planning experiments: understanding of correct ap- 
proximations and their region of applicability with the help of a consistent use 
of dimensional analysis can simplify many tasks a lot, including practical and 
non-trivial ones. Moreover, at the time being, no code capable to deal with 
transverse coherence problems has been developed at all. 



It should be noted that some theoretical attempt to follow this path has been 
proposed in j4|. Among the results of that paper is the fact that van Cittert- 
Zernike theorem could not be applied unless the electron beam divergence is 
much smaller than the diffraction angle, which is never verified in practice in 
the horizontal plane. We will show that this conclusion is incorrect. 



We organize our work as follows. After this Introduction, in Section 2 we 
present a second-order theory of coherence for fields generated by Synchrotron 
Radiation sources. In Section 3 we give a derivation of the cross-spectral den- 
sity for undulator-based sources in reduced units. Subsequently, we analyze 
the evolution of the cross-spectral density function through the beamline in 
the limit for C = 0. A particular case of quasi- homogeneous sources and its ap- 
plicability region is treated under several simplifying assumptions in Section 4. 
Effects of the vertical emittance on the cross-spectral density are discussed in 
detail in Section 5, while a treatment of some non quasi- homogeneous source is 
given in the following Section 6. Obtained results include approximate design 
formula capable of describing in very simple terms the evolution of the coher- 
ence length along the beamline in many situation of practical interest. A good 
physical insight is useful to identify possible applications of given phenomena. 
In particular in Section 7 we selected one practical application to exploit the 
power of our approach. We show that, by means of a simple vertical slit, it 
is possible to manipulate transverse coherence properties of an X-ray beam 
to obtain a convenient coherent spot-size on the sample. This invention was 
devised almost entirely on the basis of theoretical ideas of rather complex and 
abstract nature which have been described in this paper. Finally, in Section 
8, we come to conclusions. 
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Fig. 1. Light intensity from an incandescent lamp driven by a constant electric 
current. A statistically stationary wave has an average that does not vary with 
time. 

2 Second-order coherence theory of fields generated by Synchrotron 
Radiation sources 



2. 1 Thermal light and Synchrotron Radiation: some concepts and definitions 



A great majority of optical sources emits thermal light. Such is the case of the 
sun and the other stars, as well as of incandescent lamps. This kind of radia- 
tion consists of a large number of independent contributions (radiating atoms) 
and is characterized by random amplitudes and phases in space and time. The 
electromagnetic fields can be then conveniently described in terms of Statisti- 
cal Optics, a branch of Physics that has been intensively developed during the 
last few decades. Today one can take advantage of a lot of existing experience 
and theoretical basis for the descriptions of fluctuating electromagnetic fields 

Consider the light emitted by a thermal source passing through a polariza- 
tion analyzer (see Fig. 1). Properties of polarized thermal light are well-known 
in Statistical Optics, and are referred to as properties of completely chaotic, 
polarized light p, . Thermal light is a statistical random process and state- 
ments about such process are probabilistic statements. Statistical processes 
are handled using the concept of statistical ensemble, drawn from Statistical 
Mechanics, and statistical averages are performed indeed, over many ensem- 
bles, or realizations, or outcomes of the statistical process under study. 
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Polarized thermal light is a very particular kind of random process in that 
it is Gaussian, stationary and ergodic. Let us discuss these characteristics in 
more detail. 

The properties of Gaussian random processes are well-known in Statistical 
Optics. For instance, the real and imaginary part of the complex amplitudes 
of the electric field from a polarized thermal source have Gaussian distribu- 
tion, while the instantaneous radiation power fluctuates in accordance with the 
negative exponential distribution. Gaussian statistics alone, guarantees that 
higher-order correlation functions can be expressed in terms of second-order 
correlation functions. Moreover, it can be shown [5] that a linearly filtered 
Gaussian process is also a Gaussian random process. As a result, the presence 
of a spectral filter (monochromator) and a spatial filter as in the system de- 
picted in Fig. 1 do not change the statistics of the signal, because they simply 
act as linear filters. 

Stationarity is a subtle concept. There are different kinds of stationarity. 
Strict-stationarity means that all ensemble averages are independent on time. 
Wide-sense stationarity means that the signal average is independent on time 
and that the second order correlation function in time depends only on the 
difference of the observation times. However, for Gaussian processes strict and 
wide-sense stationarity coincide 0, As a consequence of the definition of 
stationarity, necessary condition for a certain process to be stationary is that 
the signal last forever. Yet, if a signal lasts much longer than its coherence 
time t c (which fixes the short-scale duration of the field fluctuations) and it is 
observed for a time much shorter than its duration cry, but much longer than 
its coherence time it can be reasonably considered as everlasting and it has a 
chance to be stationary as well, as in the case of thermal light. 

Ergodicity is a stronger requirement than stationarity. Qualitatively, we may 
state that if, for a given random process all ensemble averages can be substi- 
tuted by time averages, the process under study is said to be ergodic: all the 
statistical properties of the process can be derived from one single realization. 
A process must be strictly stationary in order to be ergodic. There exist sta- 
tionary processes which are not ergodic. One may consider, for instance, the 
random constant process: this is trivially strictly stationary, but not ergodic 
because a single (constant) realization of the process does not allow one to 
characterize the process from a statistical viewpoint. However, this is a patho- 
logic case when both the coherence time r c and the duration time of the signal 
<7t are infinite. On the contrary, a stationary process like the radiation from 
an incandescent lamp driven by a constant current has, virtually, infinite du- 
ration. In this case different ensembles are simply different observations, for 
given time intervals, of the same, statistically identical phenomenon: then, the 
concept of ensemble average and time average are equivalent and the process 
is also ergodic. 
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Fig. 2. The intensity of an X-ray beam from a Synchrotron Radiation source. 
A statistically non-stationary wave has a time-varying intensity averaged over an 
ensemble of bunches. 

Statistical Optics was developed starting with signals characterized by Gaus- 
sian statistics, stationarity and ergodicity. Let us consider any Synchrotron 
radiation source. Like thermal light, also Synchrotron Radiation is a random 
process. In fact, relativistic electrons in a storage ring emit Synchrotron Ra- 
diation passing through bending magnets or undulators. The electron beam 
shot noise causes fluctuations of the beam density which are random in time 
and space from bunch to bunch. As a result, the radiation produced has ran- 
dom amplitudes and phases. As already declared in the Introduction we will 
demonstrate that the SR field obeys Gaussian statistics. In contrast with ther- 
mal light though, Synchrotron Radiation is intrinsically non-stationary (and, 
therefore, non-ergodic) because even if its short pulse duration cannot be 
resolved by detectors working in the time domain, it can nonetheless be re- 
solved by detectors working in frequency domain. For this reason, in what 
follows the averaging brackets (...) will always indicate the ensemble aver- 
age over bunches. In spite of differences with respect to the simpler case of 
thermal light, as we will see in this paper, also Synchrotron Radiation fields 
can be described in terms of Statistical Optics. Fig. 2 shows the geometry of 
the experiment under consideration. The problem is to describe the statistical 
properties of Synchrotron Radiation at the detector installed after the spatial 
and spectral filters. Radiation at the detector consists of a carrier modulation 
of frequency u subjected to random amplitude and phase modulation. The 
Fourier decomposition of the radiation contains frequencies spread about the 
monochromator bandwidth Au m : it is not possible, in practice, to resolve the 
oscillations of the radiation fields which occur at the frequency of the carrier 
modulation. It is therefore appropriate, for comparison with experimental re- 
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Fig. 3. Reciprocal width relations of Fourier transform pairs. 

suits, to average the theoretical results over a cycle of oscillations of the carrier 
modulation. 

Fig. 3 gives a qualitative illustration of the type of fluctuations that occur in 
cycle-averaged Synchrotron Radiation beam intensity. Within some character- 
istic time, a given random function appears to be smooth, but when observed 
at larger scales the same random function exhibits "rough" variations. The 
time scale of random fluctuations is the coherence time r c . When r c <C a T 
the radiation beyond the monochromator is partially coherent. This case is 
shown in Fig. 3: there, we can estimate r c ~ Au;" 1 . If the radiation beyond 
the monochromator is partially coherent, a spiky spectrum is to be expected. 
The nature of the spikes is easily described in terms of Fourier transform the- 
ory. We can expect that the typical width of the spectrum envelope should be 
of order of Auj/uj ~ (t c u)~ 1 . Also, the spectrum of the radiation from a bunch 
with typical duration ot at the source plane should contain spikes with char- 
acteristic width Au/u ~ (c<j(T T ) _1 , as a consequence of the reciprocal width 
relations of Fourier transform pairs (see, again, Fig. 3). 



2.2 Second- order correlations in space-frequency domain 



We start our discussion in the most generic way possible, considering a fixed 
polarization component of the Fourier transform at frequency oo of the electric 
field produced at location (z ,f± ), in some cartesian coordinate system, by a 
given collection of sources. We will denote it with E±(z , r± , lo) and it will be 
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linked to the time domain field E±(z ,r± ,t) through the Fourier transform 

oo 

E±(u>) = J dtE ± (t)e iut , (1) 

— oo 

so that 

oo 

E ± (t) = — J duE^e-™ . (2) 

— oo 

This very general collection of sources includes the case of an ultra relativistic 
electron beam going through a certain magnetic system and in particular an 
undulator, which is our case of interest. In this case z Q is simply the observa- 
tion distance along the optical axis of the system and r*j_ are the transverse 
coordinates of the observer on the observation plane. The contribution of the 
fc-th electron to the field Fourier transform at the observation point depends on 
the transverse offset (l x kJyk) an d deflection angles (r] xk ,r]y k ) that the electron 
has at the entrance of the system with respect to the optical axis. Moreover, 
an arrival time t k at the system entrance has the effect of multiplying the field 
Fourier transform by a phase factor exp (iutk) (that is, in time domain the 
electric field is retarded by a time tk)- At this point we do not need to spec- 
ify explicitly the dependence on offset and deflection. The total field Fourier 
transform can be written as 

N 

E_l(z , f± , uj) = J2 E s± (rf k , Z fc , z , f± , u) exp (iut k ) , (3) 
k=i 



where ff k ,l k and t k are random variables and N is the number of electrons 
in the beam. It follows from Eq. (3) that the Fourier transform of the Syn- 
chrotron Radiation pulse at a fixed frequency and a fixed point in space is a 
sum of a great many independent contributions, one for each electron, of the 
form E s ±(ffk, h, z o, r± , uj) exp {iuit k ). For simplicity we make three assump- 
tions about the statistical properties of elementary phasors composing the 
sum, which are generally satisfied in Synchrotron Radiation problems of in- 
terest. 

1) We assume that for a beam circulating in a storage ring random variables 

— * 

t n are independent from fj n and l n . This is always verified, because the random 
arrival times of electrons, due to shot noise, do not depend on the electrons 
offset and deflection with respect to the z-direction. Eq. (3) states that the 
fc-th elementary contribution to the total E± can be written as a product of 
the complex phasors exp (iuit k ), and E s ± that, in its turn, can be written as 
a product of modulus and phase as E s ±(f) k ,lk, z o,r±o,u) =\ E s ± k \ exp(i^ fc ). 
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Under the assumption of statistical independence of t n from fj n and l n the 
complex phasors exp (iutk), and E s ± are statistically independent of each other 
and of all the other elementary phasors for different values of k. The ensemble 
average of a given function / of random variables fj n , l n and t n is by definition: 



oo 

dtr 



(f (Vni Int tn}} J dlf] xn J df]y n J dl xn j dly n J 

— oo — oo — oo — oo — oo 

xf(rf n ,l n ,t n }P(rf n ,l n ,t n } ,(4) 
where P(ff n: l n ,t n } is the probability density distribution in the joint random 

— * — * 

variables ff n , l n , t n . Independence of t n from fj n and l n allows us to write 

Pyint^nitn) = Fr] x ,l x ylxm ^xri}Fri y ,l y ijlyni lynj^titn} i (5) 



where we also assumed that the distribution in the horizontal and vertical 
planes are not correlated. Since electrons arrival times are completely uncor- 
rected from transverse coordinates and offsets, the shapes of F Vx j x , F Vy j y and 
F t are the same for all electrons. 

2) We assume that the random variables | E s±k | (at fixed frequency a;), are 
identically distributed for all values of k, with a finite mean (| E s±k |) and a 
finite second moment (| E s j_k | 2 ). This is always the case in practice because 
electrons are identical particles. 

3) We assume that the electron bunch duration <j T is large enough so that 
ujot 1: under this assumption the phases utk can be regarded as uniformly 
distributed on the interval (0, 2tt). The assumption ua T ^> 1 is justified by 
the fact that oo is the undulator resonant frequency, which is high enough to 
guarantee that ojot 1 for any practical choice of ot- 

The formal summation of phasors with random lengths and phases is illus- 
trated in Fig. 4. Under the three previously discussed assumptions we can use 
the central limit theorem to conclude that the real and the imaginary part 
of E± are distributed in accordance to a Gaussian law. Detailed proof of this 
fact is given in Appendix A. As a result, Synchrotron Radiation is a Gaussian 
random process and second-order field correlation function is all we need in 
order to specify the field statistical properties. In fact, as already remarked, 
higher-order correlation functions can be expressed in terms of second-order 
correlation functions. 

In Synchrotron Radiation experiments with third generation light sources de- 
tectors are limited to about 100 ps time resolution and are by no means able 
to resolve a single X-ray pulse in time domain: they work, instead, by count- 
ing the number of photons at a certain frequency over an integration time 
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Fig. 4. Amplitude and phase of the resultant vector (total complex amplitude) 
formed by a large number of complex phasors having random length and random 
phase. 

longer than the pulse. Therefore, for Synchrotron Radiation related issues the 
frequency domain is much more natural a choice than the time domain, and 
we will deal with signals in the frequency domain throughout this paper. The 
knowledge of the second-order field correlation function in frequency domain 

F w (z , f ±ol ,f ±o2 , u, u') = (E±(z , r ±ol ,u)E* ± (z , f ±o2 , u/)) , (6) 

is all we need to completely characterize the signal from a statistical view- 
point. For the sake of completeness it is nonetheless interesting to remark 
that it is possible (and often done, in Statistical Optics) to give equivalent 
descriptions of the process in time domain as well. First, note that the time 
domain process E±(t) is linked to E±(u) by Fourier transform, and that a lin- 
early filtered Gaussian process is also a Gaussian process (see 3.6.2). As a 
result, E±(t) is a Gaussian process as well. Second, the operation of ensemble 
average is linear with respect to Fourier transform integration. This guaran- 
tees, that the knowledge of in frequency domain is completely equivalent to 
the knowledge of the second-order correlation function between E±(z , r± a i, ti) 
and E±(z a , r± o2 , t 2 ). The latter is usually known as mutual coherence function 
and was first introduced in |9||: 



For the rest of this paper we will abandon almost entirely any reference to the 
time domain and work consistently in frequency domain with the help of Eq. 
(6) because, as has already been said, this is a natural choice for Synchrotron 



r t (z ,r_ Lol ,r_ Lo 2,ti,t 2 ) = (E_ L (z ,f_ Lo i,t 1 )E* ± (z ,f± o2 ,t 2 )) 



(7) 
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Radiation applications. In particular, as has already been anticipated, under 
non-restrictive assumptions on characteristic bandwidths of the process, it is 
possible to break the correlation function r u (u, u>') in space-frequency domain 
in the product of two factors, the spectral correlation function F w {u — u>') , 
and the cross-spectral density of the process G UJ (z , fj_ i, r± o2 , a>) jfj. The cross- 
spectral density can be studied independently at any given frequency giving 
information on the spatial correlation of the field. Subsequently, the knowledge 
of the spectral correlation function brings back the full expression for r w . 

Substituting Eq. (3) in Eq. (6) one has 



With the help of Eq. (4) and Eq. (5), the ensemble average (exp {iut k )) t can 
be written as the Fourier transform of the bunch longitudinal profile function 
F t (t k ), that is 




Expanding Eq. (8) one has 




xE* ± (ff m , l m , z , rui, u') exp - u')t m ] 




oo 




(10) 



— oo 



Using Eq. (10), Eq. (9) can be written as 



x E*j_(rj m , l m , 



m=l 
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x( E s± (rj m , I 




E*±(Vn,ln,Z ,f ±o2 ,Uj') 




where F*{uj') = F w (—uj') because F t is a real function. When the radiation 
wavelengths of interest are much shorter than the bunch length we can safely 
neglect the second term on the right hand side of Eq. (11) since the form fac- 
tor product F w {u)F u {—uj') goes rapidly to zero for frequencies larger than the 
characteristic frequency associated with the bunch length: think for instance, 
at a centimeter long bunch compared with radiation in the Angstrom wave- 
length range. It should be noted, however, that when the radiation wavelength 
of interested is longer than the bunch length the second term in Eq. (11) is 
dominant with respect to the first, because it scales with the number of par- 
ticles squared: in this case, analysis of the second term leads to a treatment 
of Coherent Synchrotron Radiation phenomena (CSR). In this paper we will 
not be concerned with CSR and we will neglect the second term in Eq. (11), 
assuming that the radiation wavelength of interest is shorter than the bunch 
length: then, it should be noted that F w {uj — uj') depends on the difference 
between u and u/, and the first term cannot be neglected. We can therefore 
write 



As one can see from Eq. (12) each electron is correlated just with itself: cross- 
correlation terms between different electrons was, in fact, included in the sec- 
ond term on the right hand side of Eq. (11), which has been dropped. It is 
important to note that if the dependence of E s ± on uj and uj' is slow enough, 
so that E s ± does not vary appreciably on the characteristic scale of F u , we 
can substitute E* ± (ff, I, z a , r± o2 , uj') with E* ± (ff,l, z ,r± o2 ,uj) in Eq. (12). The 
situation is depicted in Fig. 5. On the one hand, the characteristic scale of 
F w is given by l/crj-, where gt is the characteristic bunch duration. On the 
other hand, the bandwidth of single particle undulator radiation at resonance 
is given by u /N w , where uj is the resonant frequency and N w is the number 
of undulator periods (of order 10 2 — 10 3 ). In the case of an electron beam the 
undulator spectrum will exhibit a longer tail, as has been shown in 0|, which 
guarantees that uj /N w is, indeed, a minimum for the radiation bandwidth, 
and is the right quantity to be compared with l/oy. As an example, for wave- 
lengths of order lA, N w ~ 10 3 and a T ~ 30 ps (see Q), u a /N w ~ 2 • 10 16 Hz 
which is much larger than 1/gt ~ 3 • 10 10 Hz. As a result we can simplify Eq. 



A' 



^uj{zo,r±oi,r ±0 2,uJ,uj') = F w (uj-uj') 



m=l 




(12) to 
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F t (cJ-cj') 

* 10" 8 for hard X-rays 

> Q 

Fig. 5. Schematic illustration of the relative frequency dependence of the spec- 
tral correlation function F w (u — uj') and of the cross-spectral density function (the 
cross-power spectrum) G u (z ,r± i,r± 2,(jj) of the Synchrotron Radiation at points 
rj_ i and r_\_ o2 at frequency uj. 

T w {z , f ±ol ,f ±o2 , uj, uj') = NF w (uj - iv')G u (z , f± i, f ±o2 , a;) (13) 
where 

G u (z , f ±ol ,f ±o2 , uj) = ( E s± (rj, I, z , f ±ol ,u)E* ± (ff, I, z , f ±o2 , uj) J . (14) 

\ / ff,i 




Eq. (13) fully characterizes the system under study from a statistical view- 
point. However, in practical situations, the observation plane is behind a 
monochromator or, equivalently, the detector itself is capable of analyzing 
the energy of the photons. The presence of a monochromator simply modifies 
the right hand side Eq. (13) for a factor T{uj)T*{uj'), where T is the monochro- 
mator transfer function: 



Tc^o, f± i, r± 2, uj, uj') = NFu{lu - u)')T(u)T*(u')G u (z , f ±ol ,f ±o2 , w)(15) 

Independently on the characteristics (and even on the presence) of the monochro- 
mator, it should be noted that both in Eq. (13) and Eq. (15), correlation in 
frequency and space are expressed by two separate factors. In particular, in 
both these equation, spatial correlation is expressed by the cross-spectral den- 
sity function G u (z a , f± i, rj_ o2 , a;). In other words, we are able to deal separately 
with spatial and spectral part of the correlation function in space-frequency 
domain with the only non-restrictive assumption that uj q /N w ^> l/cr T . From 
now on we will be concerned with the calculation of the correlation function 
G U! (z ,f± i,f± 2,uj), independently on the shape of the remaining factors on 
the right hand side of Eq. (15) which can have a simple or a complicated 
structure, accounting for the characteristics of the monochromator. 
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Before proceeding with the analysis of G u (z , f± i, rj_ o2 , uj) though, let us spend 
some words on these remaining factors; the presence of a monochromator in- 
troduces another bandwidth of interest. If we indicate the bandwidth of the 
monochromator with Acj m and the central frequency of interest at which the 
monochromator is tuned with uj (typically, the undulator resonant frequency), 
then T is peaked around uj and goes rapidly to zero as we move out of the 
range (uj — Au; m /2, uj + Acj m /2). Now, if the characteristic bandwidth of the 
monochromator, Auj m , is large enough so that T does not vary appreciably on 
the characteristic scale of F u , i.e. Auj m 3> 1/or, then F^uj — uj') is peaked at 
uj = uj' . In this case the process resembles more and more a stationary process, 
although it will be still intrinsically non-stationary. Consider a signal observed 
for a time much shorter than its duration, but much longer than its coherence 
time, and such that the stationary model applies to it. Now imagine that we 
extend the observation time to a duration which is still much shorter than 
the signal duration, but long enough that we need to account for the intrin- 
sical non-stationarity of the process due to finite signal duration. In this case 
the stationary model does not apply anymore strictly. To describe this situa- 
tion, we can define a property weaker than stationarity, but nonetheless very 
interesting from a physical standpoint: quasi- stationarity. The time domain 
correlation function (that is, the mutual coherence function) can be written 
as 



oo oo 

r^f^fWx,^) = — — J dw f du'F u (u-u')T(u)T*(u') 

— oo — oo 

xG^Zo, f± i, f ±o2 , uj) exp (-iuti) exp (iuj't 2 ) . (16) 

When Au m ^> 1/or, and with the help of new variables Auj = uj—uj' and uj, we 
can simplify Eq. (16) accounting for the fact that F^uj — uj') is strongly peaked 
around Auj = 0. In fact we can consider T{uj)T*{uj')G w {z , r_j_ ol , r_j_ o2 , uj) ~ 
\T(uj)\ 2 G u} (z , f± i, r_j_ o2 , uj), so that we can integrate separately in Auj and uj 
to obtain 



00 

N r 

^t(z ,r±oi,f±o2,ti,t 2 ) = ttt^o j dAuj F^Auj) exp (-iAujt 2 ) 

—00 

00 

x J duj \T(uj)\ 2 G u) (z ,f ±ol ,f ±o2 ,uj)exp[-iuj(t 1 -t 2 )] 



= F t (t 2 ) G t (z , f ±ol , r Lo2 , t x - t 2 ) . (17) 

In other words, in the quasi-stationary case, T t {z , r± ol , f± o2 , t\, t 2 ) is split on 
the product of two factors, a "reduced mutual coherence function", that is 
Gt(z ,r±oi,f± o2 ,ti — 1 2 ), and an intensity profile, that is F(t 2 ). 
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If we now assume Au m N w /u <C 1 (that is usually true), G u (z , r± i, r± o2 , oj) = 
G u (z 01 r± i, r± o2 , u ) is a constant function of frequency within the monochro- 
mator line. In this case, contains all the information about spatial corre- 
lations between different point and is, in fact, the quantity of central interest 
in our study, but it is independent on the frequency uj. As a result we have 

r t (2 OJ r ±ol , f ±o2 , £1, t 2 ) = Ngtfa - t 2 )F t (t 2 ) G w (z a , f ±ol , f ±o2 , u Q ) , (18) 



which means that the mutual coherence function T t (z ,r± i,rj_ o2 ,ti,t 2 ) is re- 
ducible, in the sense that it can be factorized as a product of two factors, 
the first Ng t (t\ — t 2 )F t (t 2 ), characterizing the temporal coherence and the 
second G^^Zo, f± i, r± o2 , u ) describing the spatial coherence of the system 2 . 
This case is of practical importance. In fact for lA radiation we typically have 
Au m /u a ~ 1CT 4 + 1CT 5 and N w ~ 10 2 10 3 , i.e. Au m N w /uj < 1. It should 
be noted that, although Eq. (18) describes the case when Au m N w /u <C 1 
and Au m ^> 1/ut , only the former assumption is important for the mutual 
coherence function to be reducible. In fact, if the former is satisfied but the 
latter is not, from Eq. (16) one would simply have: 

Ft(z ,rxoi,rxo2,h,h) = Ngt(ti,t 2 )Gu{z ,r± i,rxo2,u ) , (19) 



that is still reducible. 

Eq. (17) contains two important facts: 

(a) The temporal correlation function, that is G t (z , r± a i, r± o2 , t\ —t 2 ), and 
the spectral density distribution of the source, that is H(z ,r± i,r± o2 ,uj) = 
N\T(u)\ 2 G u (z Q , r*j_ ol , r*j_ o2 , u>), form a Fourier pair. 

(b) The intensity distribution of the radiation pulse F t (t 2 ) and the spectral 
correlation function F^Au) form a Fourier pair. 

The statement (a) can be regarded as an analogue, for quasi-stationary sources, 
of the well-known Wiener-Khinchin theorem, which applies to stationary sources 
and states that the temporal correlation function and the spectral density are 

2 It is interesting, for the sake of completeness, to discuss the relation between G w 
and the mutual intensity function as usually defined in textbooks 0, 0| in quasi- 
monochromatic conditions. The assumption Auj m 3> 1/ct in the limit o~t — ► oo 
describes a stationary process. Now letting Auj m — ► slowly enough so that 
Aw m 3> 1/o~t, Eq- (17) remains valid while both F u and |T(u;)| 2 become approx- 
imated better and better by Dirac 5- functions, S(Auj) and 5{oj — u Q ), respectively. 
Then T t ~ G u exp[— iu> (ti — t 2 )]. Aside for an unessential factor, depending on the 
normalization of F w and |T(u;)| 2 , this relation between Tt and G^ allows identifica- 
tion of the mutual intensity function with G^ as in [1, 0] . 
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a Fourier pair. Since there is symmetry between time and frequency domains, 
a " anti" Wiener-Khinchin theorem must also hold, and can be obtained by the 
usual Wiener-Khinchin theorem by exchanging frequencies and times. This is 
simply constituted by the statement (b). 

The assumption of quasi- st at ionarity is not vital for the following of this work, 
since the cross-spectral density can be studied in any case as a function of fre- 
quency. In this respect it should be noted that, although in the large majority 
of the cases monochromator characteristics are not good enough to allow res- 
olution of the non-stationary process, there are cases when it is not allowed 
to treat the process as if it were quasi-stationary. For instance, in a partic- 
ular monochromator is described with a relative resolution of 10 -8 at wave- 
lengths of about lA, or lo ~ 2 • 10 19 Hz. Let us consider, as in ;§(], the case 
of radiation pulses of 32 ps duration. Under the already accepted assumption 
l/cr r <C u /N w , we can identify the radiation pulse duration with ax- Then 
we have Au m ~ 2 • 10 11 Hz which is of order of 2k /<jt ~ 2 • 10 11 Hz: this means 
that the monochromator has the capability of resolving the non-stationary 
processes in the frequency domain. On the contrary, also in this case, the 
mutual coherence function is reducible, in the sense specified before, because 
Auj u Juj ~ 1CT 8 and N w ~ 10 2 10 3 , i.e. Auj m N w /uj <C 1. However, such 
accuracy level is not usual in Synchrotron Radiation experiments. 

To sum up, condition u /N w l/or alone allows separate treatment of 
transverse coherence properties at a given frequency through the function 
Gu(z , f±oi, r± 2, oj). The condition for the monochromator bandwidth Au m ^> 
1/<tt defines a quasi-stationary process. If the monochromator bandwidth is 
such that Au m > lu /N w , the ratio N w / {u) &t) gives us a direct measure of the 
accuracy of the stationary approximation which does not depend, in this case, 
on the presence of the monochromator. When a monochromator is present, 
with a bandwidth Au m < u /N w , it is the ratio l/(Au m cr T ) which gives such 
a measure. The condition Auj m <C u /N w ensures, instead, that the mutual co- 
herence function of the signal is reducible, in the sense specified by Eq. (19). In 
the following we will need only the first of these conditions, u a /N w ^> 1/ctt- In 
fact, this is all we need in order to separate transverse and temporal coherence 
effects, as shown in Eq. (15). 

Once transverse and temporal coherence effects are separated one can focus on 
the study of transverse coherence through the function G^{z , f± x, r± o2 , lu). 

There exists an important class of sources, called quasi- homogeneous. As we 
will see, quasi- homogeneity is the spatial analogue of quasi-stationarity. 

In general, quasi-homogeneous sources are a particular class of Schell's model 
sources. Schell's model sources are defined by the condition that their cross- 
spectral density at the source plane (that is for a particular value of z a ) is of 
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Fig. 6. Illustrating the concept of quasi-homogeneous source. Spectral density varies 
so slowly with the position that it is approximatively constant over distances of the 
order of the correlation length A across the source. 

the form 



GJf± i,f± o2 ,uj) 



E(f± ol ,u) 



1/2 



E(r ±o2 ,u) 
xg{r±o2 - r±oi,w) , (20) 



1/2 



where g(r± o2 — r± ol , uj) is the spectral degree of coherence (that is normalized 
to unity by definition, i.e. g(0,u) = l) 3 . Equivalently one may simply define 
Schell's model sources using the condition that the spectral degree of coherence 
depends on the positions across the source only through the difference fj_ 2 — 
f± i (see jfi], 5.2.2) from which Eq. (20) follows. 

Quasi-homogeneous sources are Schell's sources obeying the following extra- 
assumption: the spectral density G w (f± ,f± ,u) at the source plane, consid- 
ered as a function of r± , varies so slowly with the position that it is approxi- 
matively constant over distances across the source, which are of the order of 
the correlation length A (that is the effective width of \G u (r± i,r± 2,w)\, see 
Fig. 6 for a qualitative illustration in one dimension). 

Since for quasi-homogeneous sources the spectral density is assumed to vary 
slowly we are allowed to make the approximation 



G w (ri i, r ±o2 , u) = I (f± i,uj) g(f ±o2 - f ±ol , uj) , (21) 

3 Sometimes, loosely speaking, we will refer to g as to "the cross-spectral density", 
or to "the field correlation function" the difference being just a normalization factor. 
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Fig. 7. Measurement of the cross-spectral density of a undulator source, (a) Young's 
double-pinhole interferometer demonstrating the coherence properties of undula- 
tor radiation. The radiation beyond the pinholes must be spectrally filtered by 
a monochromator or detector (not shown in figure), (b) The fringe visibility of 
the resultant interference pattern is equal to the absolute value of the normalized 
cross-spectral density V = \g(r± 2 ~ ^J_ol) w )|- 

where 

I(f± ol ,uj) = (\E(f ±oU Lu)\ 2 ) (22) 




fringes 




is the field intensity distribution. A schematic illustration of the measurement 
of the cross-spectral density of a undulator source is given in Fig. 7. 

We reported definitions and differences between Schell's sources and quasi- 
homogeneous sources as treated in [6j in order to review some conventional 
language. However, in our paper we will study and classify sources with the 
help of parameters from dimensional analysis of the problem. In particular, in 
the following we will introduce quantities that model, in dimensionless units, 
the electron beam dimension and divergence in the x and y directions (N x ,y 
and D x ,y, respectively). When some of these parameters are much larger or 
much smaller than unity, in certain particular combinations discussed in the 
following part of this work, we will be able to point out simplifications of an- 
alytical expressions. These simplifications do not depend, in general, on the 
fact that the source is quasi-homogeneous or not, but simply on the fact that 
some of the above parameters are large or small. It will be possible to de- 
scribe some parameter combination in terms of Schell's or quasi-homogeneous 
model, but this will not be the case, in general. In order to link the physical 
properties of certain kind of sources to a certain range of parameters we will 
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find it convenient to extend the concept of quasi-homogeneity to the concept 
of "weak quasi-homogeneity". It is better to familiarize with our new defini- 
tion already here: a given wavefront at fixed position z will be said to be 
weakly quasi-homogeneous, by definition, when the modulus of the spectral 
degree of coherence \g\ depends on the position across the source only through 
the difference r± o2 — f± i, i.e. \g\ = w(r± o2 — r±oi,u). This is equivalent to 
generalize Eq. (21) to 

\G U} (f ±ol ,f ±o2 ,uj)\ = I (r ±ol ,u})w(r ±o2 -r ±ol ,uj) . (23) 

As a remark to both the definitions of quasi-homogeneity and weak quasi- 
homogeneity, it should be noted that they only involve conditions on the cross- 
spectral densities through Eq. (21) or Eq. (23): one can apply these definitions 
to any wavefront at any position z . This is consistent with what has been re- 
marked in the Introduction: the choice of a source plane down the beamline 
is just a convention. Therefore, our definition of weak quasi-homogeneity is 
completely separated from the concept of source plane. It may seem, at first 
glance, that the definition of "weak quasi-homogeneity" is somehow artifi- 
cial but it is, on the contrary, very convenient from a practical viewpoint. 
In fact, in any coherent experiment, the specimen is illuminated by coherent 
light from some kind of aperture, or diaphragm. Think, for instance, to the 
usual process of selection of transversely coherent light through a spatial fil- 
ter, where a diaphragm is placed downstream a pinhole. Physically, when the 
modulus of the spectral degree of coherence depends only on the coordinate 
difference f± o2 — f± a i the coherence properties of the beam do not depend on 
the position of the diaphragm with respect to the transverse coordinate of 
the center of the pinhole, which we may imagine on the z axis: for instance, 
the coherence length will depend only on fj_ o2 — fj_ ol and not on the average 
position (r± o2 + rj_ ol )/2. Mathematically, as has already been said, the defini- 
tion of weak quasi-homogeneity is linked with particular combination of small 
and large parameters which will lead to simplifications of equations and to 
analytical treatment of several interesting cases. 

After having introduced the definition of "weak quasi-homogeneity" we should 
go back to the concept of usual quasi-homogeneity to describe a vary particu- 
lar feature of it: quasi-homogeneity can be regarded as the spatial equivalent 
of quasi-stationarity, as anticipated before. Exactly as the time domain has a 
reciprocal description in terms of frequency, the space domain has a reciprocal 
description in terms of transverse (two-dimensional) wave vectors. However, 
since the frequency is fixed, the ratio between the horizontal or vertical com- 
ponent of the wave vector and the longitudinal wave number is representative 
of the propagation angle of a plane wave at fixed frequency. Therefore any 
given signal on a two-dimensional plane can be represented in terms of su- 
perposition of plane waves with the same frequency and different angles of 
propagation, which goes under the name of angular spectrum. This defines an 
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angular domain which is the reciprocal of the space domain. Intuitively the 
angular spectrum representation constitutes a picture of the effects of propa- 
gation in the far zone whereas the near zone is described by the space domain. 
In this Section, an analogous of the Wiener-Khinchin theorem and its recip- 
rocal form for quasi-stationary processes has been discussed. Substitution of 
times with position vectors and frequencies with angular vectors allows to de- 
rive similar statements for the near (space domain) and far (angular domain) 
zone 0. If we use 9 = r^ jz (and #i j2 = r± ol:2 /z ) as variables to describe 
radiant intensity and cross-spectral density in the far field, and if we identify 
points on the source plane with r± (and rj_i j2 ), we have: 

(a') The cross-spectral density of the field at the source plane g(r± 2 — r±i) 
and the angular distribution of the radiant intensity 1(6) are a Fourier Pair. 

(b') The cross-spectral density of the far field g(0 2 —0 1 ) and the source-intensity 
distribution I(r±) are, apart for a simple geometrical phase factor, a Fourier 
Pair. 

The statement (b') can be regarded as an analogue, for quasi-homogeneous 
sources, of the far-zone form of the van Cittert-Zernike theorem. The statement 
(a') instead, is due to the symmetry between space and angle domains, and 
can be seen as an "anti" VCZ theorem. This discussion underlines the link 
between the VCZ theorem and the Wiener-Khinchin theorem. 

Many undulator radiation sources are quasi- homogeneous sources in the usual 
way. In the case of a quasi-homogeneous source of typical linear dimension d, 
the angular spectrum at distance z Q ^> duA/c is expected to exhibit speckles 
with typical linear dimension z (du> / c)" 1 , as illustrated in Fig. 8, which shows 
an intuitive picture of the propagation of transverse coherence. 

They generate fields which are relatively simple to analyze mathematically and 
still rich in physical features. Although both VCZ and "anti" VCZ theorem 
are based on the assumption of usual quasi-homogeneity we will see that these 
are often applicable, at least in some sense, also in the case of "weakly quasi- 
homogeneous" wavefronts. 



3 Evolution of the cross-spectral density function through the 
undulator beamline 

In our work we presented an expression for the reduced field E±_(u) = 
E±(u) exp (—iuoZo/c) of a single particle with offset I and deflection ff with 
respect to the optical axis z in an undulator. In order to derive our result, 
we used a Green's function approach to solve the paraxial Maxwell equations 
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Fig. 8. Geometry for propagation of transverse coherence in the case of a 
quasi-homogeneous source. 



for the Fourier transform of the electric field and we took advantage of a 
consistent use of the resonance approximation. The field E± differs from E± 
for a phase factor which depends on the variable z Q and on the frequency uj 
only: therefore, the use of one expression instead of the other in the equation 
for G does not change the result. In [sj], we presented results in normalized 
units in the far field zone for a particle with offset and deflection. Based on 
that work we can calculate the field in normalized units for a particle with 
offset and deflection at any distance from the exit of the undulator, where the 
center of the undulator is taken at z = 0, as specified in Fig. 9: 



1/2 

E s ± = z Q dz! 



-1/2 



z - z 



7 CX P ' 



C + 



z' + 



Wo - i0 



.(24) 



Eq. (24) is valid for the system tuned at resonance with the fundamental har- 
monic uj . This means that we are considering a large number of undulator 
periods N w 3> 1 and that we are looking at frequencies near the fundamen- 
tal and at angles within the main lobe of the directivity diagram. In this 
situation one can neglect the vertical ^/-polarization of the field with an accu- 
racy (47rA^) _1 . This constitutes a great simplification of the problem since, 
at any position of the observer, we may consider the electric field Fourier 
transform, E s ±, as a complex scalar quantity corresponding to the surviving 
x-polarization component of the original vector quantity. Normalized units 
were defined as 
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Fig. 9. Illustration of the undulator geometry and of the observation plane. 
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K being the deflection parameter, L w being the undulator length, 
/ K 2 \ ( K 2 \ 



Jo 



UJn 



4 + 2K 2 y 

Aixc~i 2 
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-Ji 



4 + 2K 2 



(26) 
(27) 



being the resonant frequency, J n the Bessel function of the first kind of order 
n, X w the undulator period, (— e) the electron charge and 7 the relativistic 

Lorentz factor. / is the normalized offset in the center of the undulator. Finally, 
the parameter C represents the normalized detuning, which accounts for small 
deviation in frequency from resonance. 
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As it is shown in Appendix B, Eq. (24) can be rewritten as 



where 



1/2 

■/ 

-1/2 



2n 



z n dz' 



-exp < 



$u + Cz' + 



ZpZ 

2(z - z>) 



(28) 



(29) 



represents the observation angle and $[/ is given by 



u 



Zo 2 



(30) 



Eq. (28) is of the form 



E s ± I C, z , 9 - j - fj | = exp 







2- 








C, z , 













(31) 



Starting from the next Section we will restrict our attention to the case C — 
for simplicity. Therefore, it may be interesting to note that in the particular 
case (7 = 0, the function S can be represented in terms of the exponential 
integral function Ei as: 



S(0,z o ,( 2 ) = exp(-^ C 2 /2)4 



iziC 



Ei 



2a2 



-l + 2i 0/ 



-Ei 



l + 2z Q/ 



(32) 



It is easy to show that the expression for the function S(-) reduces to a sinc(-) 
function as z Q ^> 1 . In fact, in this limiting case, the expression for the electric 
field from a single particle, given in Eq. (24) is simplified to 



1/2 

E s ± = exp (i&u) J dz exp { iz 

-1/2 



(33) 



Eq. (33) can be integrated analytically giving 
E s± = exp (i$u) sine f — + — J , 



(34) 
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Fig. 10. Comparison between /(£) = sinc(£ 2 /4) (solid line), the real (dashed line) 
and the imaginary (dash-dotted line) parts of /(C) = S(0,z o ,( 2 ) at z Q = 1. 

where 



Zo 



(35) 



A comparison between sinc(C 2 /4) and the real and imaginary parts of 5(0, z a , ( 2 ) 
for z = 1 is given in Fig. 10. 

Let us now go back to the general case for z Q ^ 1/2 and use Eq. (28) to 
calculate the cross-spectral density. The cross-spectral density G u is given Eq. 
(14) in dimensional units and as a function of dimensional variables. Since the 
field in Eq. (28) is given in normalized units and as a function of normalized 

— # 

variables z Q , 9 XtV and C, it is convenient to introduce a version of defined 
by means of the field in normalized units: 



G(z , 6 l ,6 2 ,C) = { E s± ICBoX-j-v] E* a ± [c,zj 2 -j-fj 



Transformation of G^ in Eq. (14) to G (and viceversa) can be easily performed 
shifting from dimensional to normalized variables and multiplying G w by an 
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inessential factor: 
G 



\KueL w Ajj 



(37) 



As a result, we can always use G in stance of G^. Substituting Eq. (28) in Eq. 
(36) we obtain: 



G(z ,9 1 ,9 2 ,C) = (S 
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x exp < i 





Expanding the exponent in the exponential factor in the right hand side of 
Eq. (38), it is easy to see that terms in l\ cancel out. Terms in 9 2 xy contribute 
for a common factor, and only linear terms in l x>y remain inside the ensemble 
average sign. Substitution of the ensemble average with integration over the 
beam distribution function leads to 



G(z ,6 1 ,6 2 ,d) = exp 



xS 




drjdl F^j ( rj, I ] exp 
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G, z , 


l---f, 




s* 




\ z " ) 







i(9 2 - 0x) • Z 

2- 

i 



C, z , | 6 2 - j- - i) 



(39) 



Here integrals dfj and in dl are to be intended as integrals over the entire plane 

— * 

spanned by the fj and I vectors. Eq. (39) is very general and can be used as a 
starting point for computer simulations. 



We already assumed that the distribution in the horizontal and vertical planes 
are not correlated, so that 7 = F~ t F~ t . If the transverse phase space is 
specified at position z Q = corresponding to the minimal values of the f3- 
functions, we can write F- , 



l x i)ySy 



F^ y F iy with 



exp 



. . 1 



2ttA, 



2D a 



2A, 
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From Eq. (25) and Eq. (29) it is easy to see that 



A 2 uL w 

U x,y — a x>,y> 



N = a 

iv z,J/ u x,y 



LU 



(41) 
(42) 



where a x>y and <7 x ',y' are the rms transverse bunch dimension and angular 
spread. Parameters N XjV will be indicated as the beam diffraction parameters 
and are, in fact, analogous to Fresnel numbers and correspond to the normal- 
ized square of the electron beam sizes, whereas D xy represent the normalized 
square of the electron beam divergences. 

Substitution of relations (40) in Eq. (39) yields, at perfect resonance ((7 = 0): 



cxp 



i (el - 1) z /2 




For notational simplicity, in Eq. (43) we have substituted the proper notation 

G(z , §i, §2, C) with the simplified dependence G(z , §i, e 2 ) because we will be 

treating the case C — only. Consistently, also S[z Q , (9 — 1/ z Q — fj) 2 ] is to be 

understood as a shortcut notation for S[C,z Q , (9 — l/z — f)) 2 ] calculated at 
C = 0. 



4 Undulator radiation as a quasi-homogeneous source 



When describing physical principles it is always important to find a model 
which provides the possibility of an analytical description without loss of es- 
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sential information about the feature of the random process. 

In order to get a feeling for some realistic magnitude of parameters we start 
noting that the geometrical emittances of the electron beam are simply given 
by t x ,y — u x,y a x' ,y' ■ Here they will be normalized as e x>y = 2ire Xty /X. Then 
a xy = Pxy e x,y^ where f3° y are the minimal values of the horizontal and verti- 
cal betatron functions. In this paper we will assume that the betatron func- 
tions will have their minimal value at the undulator center. Therefore we have 
e x,y = a l', y 'Px,y or > in normalized units , = D X: yP Xt y, where /3 XjJ/ = $° X JL W . 
Equivalently we can write e x>y = \J D x ^N x ^ y . It follows that N Xjy = e Xty j3 X;y . 
Now taking A = 1 A, e x = 1 + 3 nm, e y = I0~ 2 e x and (3° = 1CT 1 + 10L W one 
obtains, in normalized units, e x = 10 2 + 3 ■ 10 2 , e y = 1 + 3 and (3 X = 10 1 + 10: 
therefore D x >> 1, N x >> 1. This is always the case in situations of practical 
interest, with N x which may range from values much smaller to much larger 
than D x . 

Assuming D x ^> 1 and N x 1, independently on the values of D y and N y , 
introduces simplifications in the expression for the cross-correlation function 
and allows further analytical investigations. As we will see, in particular, a 
model of the electron beam based on these assumptions contains (but it is not 
limited to) the class of quasi-homogeneous sources discussed in Section 2.2. 

In the next Sections 4.1 and 4.2 and later in Section 5, we will see what are the 
conditions in terms of the dimensionless parameters N XjV and D x y for some 
undulator radiation wavefront at position z a , to be quasi-homogeneous in the 
usual and in the weak sense (according to the definition in Section 2.2), we 
will justify the introduction of the concept of weak quasi- homogeneity itself 
and we will discuss the applicability regions of the VCZ (and "anti" VCZ) 
theorem. Then, in Section 6, we will also discuss some case characterized by 
non weakly quasi-homogenous fields. 

4-1 A simple model 

To provide a first analysis of the problem we adopt some simplifying assump- 
tions that are only occasionally met in practice. 

As already assumed vertical emittance is much smaller than horizontal emit- 
tance. For notational simplicity we will make the assumptions N y 1 and 
D y 1. This means that we theoretically assume fj y <^ 1 and 1^ < 1. As a 
result, the terms in fj y and l y can be neglected in the S(-) term on the right 
hand side of Eq. (43). Although this model includes obvious schematization 
it is still close to reality in many situations, and it is only to be considered 
as a provisory model for physical understanding to be followed, below, by 
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more comprehensive generalizations. In this Section we will restrict our atten- 
tion at the correlation function for y \ = 9 y2 that is on any horizontal plane. 
Here again, for notational simplicity, we will substitute the proper notation 
G(z a , 9 x i, 0, 9 x2 , 0) with G(z Q , 9 xl , 9 x2 ). Eq. (43) can be greatly simplified lead- 
ing to 



G(z , 9 x i, 9 x2 ) 



1 



oo 

x / //,. exp 



oo 

J dfj x 



-oo 



2D, 



dl x exp 



cxp 



2N X 



cxp 



- 2 x2 ) z /2 
i(9 X 2 — 9 x i)l x 



xlS 



o, (dxi — lx/ z — VxJ S* z , (§ x2 — l x j z Q — rj x ^J 



.(44) 



Let us now introduce 
9 x i — 9 x2 



A9 



9 x i + 9 x2 



(45) 
(46) 



With this variables redefinition we obtain 



ex. 

/ 



x / dl x exp 



G(z ,9, A9) = 

< ii \ 



2N X 



; exp (i29A9zn) / dfi x exp ( 

S z , (0 + AO- l x j z - rjxj 



cxp 



-2iA9L 



xS* 



z ,{9-A9-l x /z -f,x) 2 ^ .(47) 



A double change of variables f) x — >■ fj + 9 followed by l x /z — > (f) — fj yields 



_ „ exp li26A9z ) °f ( (f, + e) 2 ~ ^ 

G{z ,9,A9) = \_ > I dfiexp{- Kl ^ ) + 2iA9z r) 



2nJDN/z 2 j 



x 



J # exp (- ^ts.fT ) exp(-2iA9z, 



xS* 



2N/zl 

z ,{4>-A9f]s\z ,{4> + A9) 2 ] .(48) 



where we have posed D = D x and N — N x for notational simplicity. Eq. (48) 
can be also written as 
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_ A exp(i26A9z ) 
G(z D , 9, A9) = v ^ = ' exp 



x 



2ir^/DN/z 2 
x exp (-2iA9z $) S 



exp 




, 2DN/S* ' ■ N/z 2 ' D 
The integral in rj can be performed analytically thus leading to 



.(49) 



exp (i29A9z^j 



x 



J 4 



exp 







x exp 
2 1 20 



^2n(N/z 2 + D 

# 2 + ANA9 2 D + Ai(N/z )9A9 

2(iV/4 2 + t>) . 

2i(N/z )A9, , 

S* \z , (0 - AO) 

Zo ,(4> + A9) 2 ] .(50) 



2(N/z 2 + D) 



xS 



It is important to remember again that an asymptotic formula for z ^> 1 can 
be obtained from Eq. (50) simply substituting S[z D , (0± A#) 2 ] with sinc[(0± 
A9) 2 /4]. Then, it is easy to understand that S is bound to go to zero for values 
of (0± A#) 2 larger than unity, exactly as the asymptotic terms in sinc(-) would 
do. In fact, once C is set to zero, S depends parametrically on the normalized 
distance z Q alone, that is S — S\z , (0± A#) 2 ], and gives the previously found 
asymptotic expression of sinc[(0 ± A9) 2 /4] in the limit for z 1. Since here 
z is supposed to be at least of order unity (z Q > 1/2), we can conclude that S 
must be different from zero only for values of (0 ± A9) 2 of order unity (to be 
more precise, for values (0 ± A9) 2 ~ 4, (0 ± A#) 2 /4 being the arguments of 
the sine function) as it can be seen, for instance, from Fig. 10 for a particular 
case. Thus Eq. (50) and its asymptotic equivalent for z ^> 1 share the same 
mathematical structure. 

Let us now introduce the non-restrictive assumptions: 



7V> 1 , £>> 1 
and define 



(51) 



(52) 



The physical interpretation of A follows from that of a/z : A is the dimen- 
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sionless square of the apparent angular size of the source at the observer point 
position, calculated as if the source was positioned at z — 0. If N 1 and 
D > 1 we have (2Az 2 D) /(A + D) ^> 1 for any value of z Q and any choice of 
N and D. As a result, from the exponential factor exp [— 2Az 2 A9 2 D/(A + D)] 
outside the integral sign in Eq. (50) we have that G u {z o ,0, AO) is different 
from zero only for AO <C 1. Then we can neglect terms in A6 in the factors 
S(-) within the integral sign thus getting 



_ A exp (i20z o A§) 
G(z , 0, AO) = j V = Jr exp 
2n(A + D) 

2 + 206 ' 



x 



2 + AAz 2 o A0 2 D + AiA0z o A0 
2(A + D) 



exp 



2{A + D) 



exp 



2(f)Az A9 
A + D 



—i- 



s 



(53) 



The maximal value of 9 is related with the width of the exponential function 
exp [— 9 2 j (2 A +2D)] outside the integral sign in Eq. (53). It follows that in the 
limit for D 1 we can neglect the exponential factor exp [— (0 2 + 2<f)9) / (2 A + 2D)] 
within the integral sign: in fact, its argument assumes values of order unity for 
(f)^> 1, but the factor | S[z ,(f) 2 ] | 2 cuts off the integrand for > 1. Therefore 
Eq. (53) can be simplified as follows: 



_ A exp(i29z A§) 
G{z a , 9, A9) = ; v = J - exp 



2it(A + D) 
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x exp 



2ADz 2 A9 2 
A + D 



-i oo 








■p exp 
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2{A + D) 
2A 



exp 



A + D 



z A9 



2iA9z A9 
A + D 



z ,(p 2 -(54) 



The integral in Eq. (54) is simply the Fourier transform of the function 
f(4>) =| S[z ,(p 2 } | 2 with respect to the variable -2Az A9/(A + D). Since 
the function f((f>) has values sensibly different from zero only as is of order 
unity or smaller, its Fourier Transform will also be suppressed for values of 
2Az \A9\/(A + D) larger than unity, by virtue of the Bandwidth Theorem. 
This means that the integral in Eq. (54) gives non-negligible contributions 
only up to some maximal value of \A9\: 



(55) 



On the other hand, the exponential factor outside the integral in Eq. (54) will 
cut off the function G around some other value 

( + ) . (56) 



34 



It is easy to see that, for any value of z„, \A0\ max |A#| max2 . In fact we have 
I Af)\ 

i-Al/ max 



\A6 



DyJl + (D/A) > > 1 , 



(57) 



max2 



in the limit for D ^> 1. As a result the Fourier transform in Eq. (54) is 
significant only for values of the variable — 2Az o A0j (A + D) near to zero and 
contributes to G only by the inessential factor 



/ 



dcj) S z Q , 4> 



constant 



(58) 



In order to use the correlation function G for calculation of coherence length 
and other statistical properties, one has to use the spectral degree of coherence 
g, which can be presented as a function of and AO instead of x± o2 and x± i. 



g (5, AO) 



G (e, AO) 



E s± (0 + AO) 



2\V2 / n 2 \ V2 • 

E[0- AO) 



(59) 



From Eq. (54) we obtain : 



g(z , 0, AO) = exp (i20z o A0) exp 



2iA0z o A0 

A + b 



exp 



2ADz 2 o A0 2 
A + D 



(60) 



In the asymptotic limit for a large value of z a , A 1, Eq. (60) can be simplified 
to 



g(z , 0, AO) = exp (i20z o A0) exp 



2iA0z o A0 

b 



exp 



-2AzlA0 2 



(61) 



From Eq. (54) it is easy to see that the region of interest where the field 

intensity is not negligible is when < \[~b. Therefore, Eq. (61) can be further 
approximated to 



g(z , 0, AO) = exp (i20z o A§) exp 



-2Az 2 o A0 2 



(62) 



It is interesting to calculate the transverse coherence length £ c as a function of 
the observation distance z . For any experiment, complete information on the 
coherence properties of light are given by the function g. When calculating the 
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coherence length one applies a certain algorithm to g thus extracting a single 
number. This number does not include all information about the coherence 
properties of light, and the algorithm applied to g is simply a convenient 
definition. Then, in order to calculate a coherence length one has, first, to 
choose a definition among all the possible convenient ones. In this paper we 
will simply follow the approach by Mandel, originally developed for the time 
domain, but trivially extensible to any domain of interest, in our case the 
angular domain. The coherence length, naturally normalized to the diffraction 
length yj ' L w c/u is defined as 

oo 

2 J \g(A6)\ 2 d(z A6) , (63) 

— oo 



where the factor 2 in front of the integral on the right hand side is due to the 
fact that we chose Mandel's approach and that our definition of A9 differs of 
a factor 1/2 from his definition. Performing the integration in Eq. (63) with 
the help of Eq. (60) yields: 

/ 1 IX 1 / 2 

Uzo) = ^ ( 3 + s ) . (64) 



4-2 Discussion 



The coherence length in Eq. (64) exhibits linear dependence on z , that is 
£ c — ► \Jn/N z while for z — ► 1/2 that is at the end of the undulator, it 
converges to a constant £ c — ► [n / (AN) +n / D} 1 / 2 . Eq. (64) and its asymptotes 
are presented in Fig. 11 and Fig. 12 for the case N = 10 3 , D = 10. It is evident 

that at the exit of the undulator, £ c ~ l/\[f), because N 3> D. On the other 

hand, horizontal dimension of the light spot is simply proportional to \fti 
as it is evident from Eq. (60). This means that the horizontal dimension of 
the light spot is determined by the electron beam size, as is intuitive, while 
the beam angular distribution is printed in the fine structures of the intensity 
function, that are of the dimension of the coherence length. In the limit for 
z Q ^> 1 the situation is reversed. The radiation field at the source can be 
presented as a superposition of plane waves, all at the same frequency u , but 
with different propagation angles with respect to the z-direction. Since the 
radiation at the exit of the undulator is partially coherent, a spiky angular 
spectrum is to be expected. The nature of the spikes is easily described in 
terms of Fourier transform theory, in perfect analogy with what has been said 
about the frequency spectrum in Section 2.1. From Fourier transform theorem 
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Fig. 11. Coherence length £ function of z Q and asymptotic behaviors for 

z — ► 1/2 and £ G > 1. Here JV = 10 3 and i) = 10. 

or, directly, from Eq. (54) or from geometrical optics arguments we can expect 
an angular spectrum envelope with Gaussian distribution and rms width of 




Also, the angular spectrum should contain spikes with characteristic width 

1/ \//v", as a consequence of the reciprocal width relations of Fourier transform 
pairs (see Fig. 13). This can be seen realized in mathematical form from the 
expression for the cross-spectral density, Eq. (60) and from the equation for 
the coherence length, Eq. (64). Since N ^> 1, the horizontal width of the 
coherence spot is much smaller than the vertical one. 

It is also important to remark that the asymptotic behavior for A <C 1 of g in 
Eq. (62) and | c 




are direct application of van Cittert-Zernike theorem. In fact, the last expo- 
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Fig. 12. Enlarged view of the initial part of Fig. 11. 

nential factor on the right hand side of Eq. (61) is simply linked with the 
Fourier transform of Fj (l x ). We derived Eq. (61) for N ^> 1 and D 3> 1, with 

z 2 i) 3> N: in non-normalized units these conditions mean that the VCZ the- 
orem is applicable when the electron beam divergence is much larger than the 
diffraction angle, i.e. a x , ^> A/ (2nL w ), the electron beam dimensions are much 
larger than the diffraction size, i.e. a x XL w /2tc, and (a X 'Z ) 2 ^> a x . On the 
contrary, authors of state that, in order for the van Cittert-Zernike theorem 
to be applicable, "the electron-beam divergence must be much smaller than 
the photon divergence", that is our diffraction angle, i.e. a x > <C \/(2tcL w ) 
(reference page 571, Eq. (57)). Our derivation shows that this conclusion 
is incorrect. 

In (if (paragraph 5.6.4) a rule of thumb is given for the applicability region 
of the generalization of the VCZ theorem to quasi- homogeneous sources. The 
rule of thumb requires z Q > 2dA/\ where d is "the maximum linear dimen- 
sion of the source" , that is the diameter of a source with uniform intensity 
and A "represents the maximum linear dimension of a coherence area of the 
source". In our case d ~ 2a x , since a x is the rms source dimension, and from 
Eq. (64) we have A = £ c ~ A/ \2y/iia x i). The rule of thumb then requires 
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Fig. 13. Physical interpretation of the generalized van Cittert-Zernike theorem. If 
the radiation beyond the source plane is partially coherent, a spiky angular spectrum 
is expected. The nature of the spikes in the angular spectrum is easily described 
in Fourier transform notations. We can expect that typical width of the angular 
spectrum should be of order (wA/c)" 1 , where A is the typical linear dimension 
of spatially random intensity fluctuations. Also an angular spectrum of the source 
having transverse size d should contain spikes with typical width of about (cod/c)^ 1 , 
a consequence of the reciprocal width relations of Fourier transform pairs. 

z Q > 2a x / (i/ncTx'): in dimensionless this reads z Q > J N/D. This is paramet- 
rically in agreement with our limiting condition z 2 D ^> N, even though these 
two conditions are obviously different when it come to actual estimations: 
our condition is, in fact, only an asymptotic one. To see how well condition 

%o ^ yN/D works in reality we might consider the plot in Fig. 11. There 
N = 10 3 and D = 10 so that, following [^| we may conclude that a good con- 
dition for the applicability of the VCZ theorem should be z > 10. However 
as it is seen from the figure, the linear asymptotic behavior is not yet a good 
approximation at z Q ~ 10. This may be ascribed to the fact that the deriva- 
tion in is not generally valid, but has been carried out for sources which 
drop to zero very rapidly outside the maximum linear dimension d and whose 
correlation function also drops rapidly to zero very rapidly outside maximum 
linear dimension A. 



However, at least parametrically, the applicability of the VCZ theorem in 
the asymptotic limit z 2 D 3> N can be also expected from the condition 
z Q > 2dA/X in 0. In other words, with the help of our approach we were 
able to specify an asymptotic region where the VCZ theorem holds. Such 
a region overlaps with predictions from Statistical Optics. Statistical Optics 
can describe propagation of the cross-spectral density only once it is known 
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at some source plane position. Our treatment allows us to specify the cross- 
spectral density at the exit of the undulator, but it should be noted that we 
do not need to use customary results of Statistical Optics and propagate the 
cross-spectral density from the exit of the undulator in order to obtain the 
cross-spectral density at some distance along the beamline. In fact our ap- 
proach, which consists in taking advantage of the system Green's function in 
paraxial approximation and, subsequently, of the resonance approximation, 
allows us to calculate the cross-spectral density directly at any distance from 
the exit of the undulator. 

Let us now consider the structure of Eq. (60) and discuss the meaning of the 
phase terms in 6 AO. These are important in relation with the condition for 
quasi-homogeneous source: their presence couples the two variables 6 and A9 
and prevents the source to be quasi-homogeneous at any given value of z Q 4 , 
unless they compensate each other in some parameter region. 

Let us discuss the limit, A ^> 1. We may consider two subcases. First, consider 
A 3> D ^> 1. In this case, inspection of Eq. (60) shows that the two phase 
terms compensate and the source is quasi-homogeneous, because the cross- 
spectral density is factorized in a function of 9 and a function of A6. It should 
be noted that if condition A 3> 1 is not satisfied at the exit of the undulator, 
where z a ~ 1, then it is never satisfied. If N ^> D ^> 1 we have a quasi- 
homogeneous source at the exit of the undulator. 

Second, consider D N 3> 1- This correspond to a situation with a low 
value of the normalized betatron function in the horizontal direction. Figure 
14 shows a numerical example with e x = 100 and f3 x = 0.3 that is iV = 30 and 
D = 300: the value for the horizontal betatron function is similar to the low-/3 
case reported at page 12, Table 2.2.2 in []]], where (3 X = 1.3 m for a 5m-long 
insertion device. The value e x = 100 corresponds to a wavelength of about 0.6A 
for the PETRA III case. When D > # > 1 no compensation of the phase 
terms in Eq. (60) is possible, not even at the exit of the undulator. In this case, 
whatever the value of z Q we can never have a quasi-homogeneous wavefront. 
This constitutes no problem. Simply, the wavefront is not- quasi-homogeneous 
in this case. However, we may interpret the situation by saying that a "virtual" 
quasi-homogeneous source placed in the center of the undulator would result in 
the non-homogeneous source described by Eq. (60) at the exit of the undulator. 
Although it physically makes no sense to discuss about Eq. (60) inside the 
undulator, the "virtual" source analogy is suggested by the fact that setting 
z D = in Eq. (60), both phase terms become zero. 



4 Note, again, that the definition of "source plane" is just conventional. One may 
define the source plane at the exit of the undulator, that is at z Q = 1/2, but there 
is no fundamental reason for such a definition: one may pick any value of z Q as the 
source position. 



40 



1,50 




Fig. 14. Coherence length £ c as a function of z Q in the case D 3> -/V S> 1. In 
particular, here, iV = 30 and D = 300. The linear dependence on z Q starts already 
from the exit of the undulator, because D 3> N ■ 



Note that in general, whatever the values of N 3> 1 and D>1, one never has 
quasi- homogeneous sources in the limit for A <C 1. In fact, in the asymptotic 
A <C 1 only the phase factor exp(i29A9z a ) contributes, which couples 9 and 
A9. Such a factor is connected with phase of the field from a single electron in 
an undulator in the far zone, oj(x 2 + y 2 )/(2cz ), which represents, in paraxial 
approximation, the phase difference between the point (x , y , z Q ) and the point 
(0, 0,z o )' in the asymptotic for large values of z , the electric field generated 
by a single electron with offset and deflection in an undulator has a spherical 
wavefront (see 0). When one calculates the field correlation function at two 
different points, he ends up with a contribution equal to the difference (due 
to complex conjugation) between u)(x 2 o2 + y 2 2 )/(^cz ) and oj(x 2 o1 +y 2 1 )/(2cz ), 
that for the vertical (and separately, the horizontal) direction gives exactly the 
shift 29A9z in normalized units. It should be noted that such a reasoning is 
not limited to Synchrotron Radiation sources, but it is quite general since, as 
already discussed, it relies on the fact that the wavefront of a single radiator 
(in our case, an electron) produces a spherical wavefront in the far field. This 
is, for instance, the case of thermal sources as well. In other words, if the far 
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field radiation of a quasi-homogeneous source is taken as a new source, that 
new source will never be quasi-homogeneous. 

A common property of all situations with N 3> 1 and D 3> 1 is that, for 
any value of z , the modulus of g, i.e. \g\, is always independent on 6. More- 
over, it is always possible to apply the VCZ theorem starting either from 
a virtual quasi-homogeneous source placed at center of the undulator when 
D >> N >> 1, or otherwise from a real one placed at the exit of the undulator 
(or at any other position z Q close enough to the exit of the undulator to guar- 
antee a quasi- homogeneous wavefront). These observations suggest to extend 
the concept of quasi-homogeneity, and introduce the new concept of "weak 
quasi-homogeneity" as discussed before. With respect to the new coordinates 
9 and A9, a given wavefront at fixed position z is said to be weakly quasi- 
homogeneous when \g\ is independent of 6. With this new definition at hand 
we can restate some of our conclusions in a slightly different language. We have 
seen that in the far field, when the VCZ theorem holds, the wavefronts are 
weakly quasi-homogeneous, but never quasi-homogeneous in the usual sense. 
In the case Af > D > 1 we pass from quasi- homogeneous wavefronts (in 
the usual sense) in the near field to weakly-quasi homogeneous wavefronts 
(but not quasi-homogenous in the usual sense) in the far field. Note that the 
wavefronts are always weakly quasi- homogeneous, even during the transition 
from near to far zone. In the case D ^> N ^> 1 instead, the VCZ theorem is 
applicable already from the exit of the undulator, as it can be seen from Fig. 
14, and the wavefront is not quasi-homogeneous in the usual sense, but still 
weakly quasi-homogeneous from the very beginning. 

The weak quasi-homogeneity of the wavefronts at any value of z Q , i.e. the fact 
that | g | is independent of 9 for any value of z Q guarantees that the plot in 
Fig. 14 is universal. It should be noted that this fact depends on the choice 
of large parameters N 3> 1 and D 3> 1, but it is also strictly related with the 
Gaussian nature of the electron distribution in angles and offsets, that is a 
well-established fact for storage-ring beams. If angles or offsets were obeying 
different distribution laws, in general, one could not perform the integral in fj in 
Eq. (48) and, in general, \g\ would have shown a dependence on 9: our notice- 
able result is linked with the properties of the exponential elementary function. 
However, it should be clear that even in the case when angles or offsets were 
obeying different distribution laws, i.e when the plot in Fig. 14 is not univer- 
sal, we could have situations when wavefronts are quasi-homogeneous in the 
usual sense near the exit of the undulator and are weakly quasi-homogeneous 
in the far field limit, but not along the transition between these two zones. A 
more detailed discussion of this issue will be given in Section 5, where we will 
be discussing conditions for the source to be quasi- homogeneous. 

Another remark to be made pertains the applicability of the VCZ theorem. As 
we deal with a quasi- homogenous source (in the usual sense) the knowledge of 
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1(9) and g(A9) in the far zone allow, respectively, the calculation at the source 
plane of g(Ax) through the "anti" VCZ theorem and of I(x) through the VCZ 
theorem (here we consider only one dimension, the horizontal one x). Viceversa 
the knowledge of / and g at the source allow calculation of I and g in the 
far field. In terms of intensity, all information regarding wavefront evolution 
(assuming a quasi-homogeneous source, in the usual sense) is included in 1(9) 
and I(x). For instance, the knowledge of 1(9) allows calculation of g(Ax) 
at the source plane through the "anti" VCZ theorem. Then, the knowledge 
of g(Ax) and I(x) allow the calculation of the cross-spectral density, which 
can be propagated at any distance, and allow to recover g(A9). So, complete 
characterization of the undulator source is given when 1(9) and I(x), when 
the source is assumed quasi-homogenous. 



Yet we have seen that, when the electron distribution in angles and offsets are 
Gaussian and N 3> 1 and D >> 1, the VCZ theorem holds also in the case D 3> 
N ^> 1, when the source is non quasi- homogeneous in the usual sense. We have 
seen that the situation can be equivalently described with the help of a virtual 
quasi-homogeneous source in the middle of the undulator. However, such an 
interpretation is only valid a posteriori. For the case N ^> 1 and D >> 1 there is 
a non quasi- homogeneous wavefront at the undulator exit; before our approach 
was presented one would have concluded that the VCZ theorem cannot be 
applied, since the spectral degree of coherence does not form a Fourier pair 
with the intensity distribution at the undulator exit. Our approach is based 
on the simplification of mathematical results through the use of small and 
large parameters and subsequent understanding and interpretation of these 
simplified results: in our analysis we were never limited to the treatment of 
quasi-homogenous cases alone. 



As a closing remark about the coherence length we like to draw the reader's 
attention on the fact that the dimensional form £ c of the coherence length, 
given in normalized units by £ c in Eq. (64), does not include the undulator 
length. This is to be expected since, in the limit N 3> 1 and D 3> 1, the 
typical size and divergence of the electron beam are much larger than the 
diffraction size \J cL w /u and angle \Jc~j (uL w ), which are intrinsic properties 
of the undulator radiation. As a result, in this limit, the evolution of the 
radiation beam is a function of the electron beam parameters only, and does 
not depend on the undulator length. In the following Section 5, where we 
will extend our model to a two-dimensional case, we will see that the quasi- 
homogeneous approximation is valid in many practical situations, but we will 
have to account for diffraction of undulator radiation in the vertical direction. 
In this case the dimensional coherence £ c length will be a function of the 
undulator length L w as well. 



43 



5 Effect of the vertical emittance on the cross-spectral density 



Up to now we were dealing with the field correlation function g within the 
framework of a one-dimensional model. 

In fact we considered the limit for e y (3 y <C 1 and e y /j3 y <C 1 and we calculated 
g for y i = 9 y2 = and C = 0, so that our attention was focused on coherent 
effects in the horizontal direction. We will now extend our considerations to a 
two-dimensional model always for (7 = 0. This can be done by a straightfor- 
ward generalization of Eq. (50) which can be obtained from Eq. (43) following 
the same steps which lead to Eq. (50), but this time without assumptions on 
N y , D y , 6 y i and 9 y2 . Finally, at the end of calculations, our final expression 
for G should be normalized to 




As has already been seen in Section 4, after normalization to XV we will obtain 
the spectral degree of coherence g. With this in mind we will neglect, step after 
step, unnecessary multiplicative factors that, in any case, would be finally 
disposed after normalization of the final result. Retaining indexes x and y in 
our notation we obtain 



G(z , 6 X , 9 y , A9 X , A9 y ) = exp [i2 (9J A9 X + 9 y z A9 y ) 



x exp 



x exp 



AA^zlA9lD^ + AiA x 9 x z A9 



9J + 4A y z^A9yD y 



2(A X + D x ) 

UA y 9 y Z A9y 



x J d(f) x J d(p y exp 



2(Ay + Dy) 

4> 2 X + 20^ (9 X + 2iA x z Q A9 x ) 



x exp 



2(A X + D X ) 

4> 2 y + 2^y (9y + 2iAyZ Q A9 l 



xS* 



2(Ay + Dy) 

z , (4 - A9 x f + {4> y - A9 y ) 2 S z , (4 + A9 x f + (0, + A9 y ) 2 } . (67) 



We will still assume N x ^> 1 and 4 » 1. This allows to factorize the right 
hand side of Eq. (67) in the product of contribution depending on horizontal 
(9 X , A9 X ) coordinates only with a second depending on vertical (9 y , A9 y ) co- 
ordinates only. In fact, from the exponential factor outside the integral sign 
in Eq. (67) it is possible to see that the maximum value of A9 2 £ is of order 
(A x + D x ) / (A x D x z 2 ) < 1. Asa result, A9 X can be neglected inside the S func- 
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tions in Eq. (67). Moreover, since D x ^> 1 one can also neglect the exponential 
factor in (f) 2 + 2<fi x x inside the integral sign. This leads to 



G(z , 9 X , B y , A9 X , A9 y ) = exp [i2 (9 x z A9 x + 9 y z A9 



x exp 



exp 



61 + AA x z 2 A9 2 x D x + 4iA x 9 x z A9 x 



x 



J d(f) 2 



exp 



i(f>x 



2A x z Q A9 x 



A x + D, 



d<f> y exp 



2(A X + D x ) 
9J + AA y z 2 A6 2 y b y + AiA y 9 y z A9 y 
2{A y + D y ) 

4> 2 y + 24>y (9y + 2iAyZ A9y 



2(Ay + Dy) 



z , 4> 2 X + (0, - A^) 2 5 z , 2 + (0, + A6 y ) 2 . (68) 



Based on the same reasoning in Section 4.1, that we repeat here for complete- 
ness, we can also neglect the phase factor in <p x under the integral in d(j) x in 
Eq. (68). Such integral in d(f) x is simply the Fourier transform of the function 



d(j) y exp 



\\ + 20y (9y + 2iAyZ A9y) 
2(i, + by) 



xS* 



Zo , 4> 2 X + (0, - ABy) 2 S z , 4> 2 X + (4 + A6 y ) 2 . (69) 



with respect to the variable — 2A x z A9 x /{A x + D x ). In the argument of the S 
functions on the right hand side of Eq. (69), 2 is always summed to positively 
defined quantities. This remark allows one to conclude that f((f) x ) has values 
sensibly different from zero only as <fi x is of order unity or smaller. Therefore, its 
Fourier Transform will also be suppressed for values of 2A x z o \A0 x \/(A x + D x ) 
larger than unity, by virtue of the Bandwidth Theorem. This means that the 
integral in d(p x in Eq. (68) gives non-negligible contributions only up to some 
maximal value of \A9 X \: 



\A9\ 



2z n 



1 + 



(70) 



On the other hand, the exponential factor outside the integral in Eq. (68) will 
cut off the function G around some other value 



IA0I 



1 



xmax2 



2z Q \D 



1 1 
+ — 



1/2 



(71) 



It is easy to see that, for any value of z Q , lA^max 3> \A9\ xmax 2- In fact we 
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have 



|^^|:rmax2 



fb x y/l + (D x /A x ) > \[t> x > 1 , 



(72) 



in the limit for D x > 1. As a result the Fourier transform in Eq. (68) is 
significant only for values of the variable — 2A x z A9 x j (A x + D x ) near to zero. 
As a result we obtain the following equation for G: 



G(z Q , 9 X , 9 y , A9 X , A9 y ) = exp [i2 [9 x z A9 x + 9 y z a A9 y 
9 2 x + AA x z 2 A9 2 x D x + AiA x 9 x z A9 x 



x exp 



exp 



2(A X + D x ) 
g + AA y z 2 A9yb y + AiA y 9 y z A9 y 



oo 

x J d(j) y exp 

— oo 



2(A y + D y ) 

+ 20y (By + 2iAyZ A9y) 



2(A y + D y ) 

x / d<t> x S* z o J 2 x + (0 y - A9 y f]s\z Al + {4> y + A9 y ) 2 } ,(73) 



where horizontal and vertical coordinates are obviously factorized. 

Eq. (73) has been derived assuming N x ^> 1 and D x ^> 1. Note that assuming 
setting A9 y = 9 y = one can obtain Eq. (60) from Eq. (73). This proves that 
Eq. (60) has a wider range of validity than that for iV^ < 1 and D y <^ 1 (as 
the reader will remember, these assumptions were made just for notational 
simplicity). In fact, as N x ^> 1 and D x ^> 1 horizontal and vertical direction 
factorize and the horizontal factor is always that in Eq. (60), independently 
on N y and D y . 

Under one of the two extra assumptions N y ^> 1 or D y 1, Eq. (73) can be 
further simplified and often describes a weakly quasi-homogeneous wavefront 
according to the definition given at the end of Section 4.1. At the end of 
the Section we will show that as D y > 1 we always have a weakly quasi- 
homogeneous wavefront, for any value of N y . It will also be seen that the same 
applies when N y ^> 1 and A y <^ 1 (far field) or A y 1 (near field) for any 
value of t) y . However, as N y ^> 1, D y < 1 and A y ~ 1 we have an intermediate 
region between the near and far region were, in general, wavefronts are not 
quasi- homogeneous, not even in the weak case. 



For simplicity of discussion we will set 9 X = 9 y 
(73) 



thus obtaining from Eq. 
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G(z ,A9 x ,A9 y ) = exp 



x 



2A x %A6lD x 
(A x + t) x ) 

oo 

J d<t> y exp 



cxp 



2A y zlAdlD y 
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+ 20 y (2iA y z A6 y ) 
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J d(f) x 
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z Q , 2 + (0, - AOyf S z , 4> 2 X + (0, + A^) 2 



(74) 



which is easier to manipulate. It should be reminded that, if one is interested in 
ascertaining the weak quasi- homogeneity of a wavefront, one has to deal with 
the full Eq. (73). Moreover, it should be noted that on the one hand, within 
the weak quasi-homogeneous case, Eq. (74) is quite general and we can extract 
from it all important information on the transverse coherence independently 
on the values of 9 X and 9 y . On the other hand though, in the case the weakly 
quasi- homogeneous assumption fails, Eq. (74), e.g. when N y >> 1, D y < 1 and 
A y ~ 1 as we will see, \g\ depends on 9 y and the study of Eq. (74) has a more 
restricted range of validity, namely for the particular value of 9 y = 0. 



In Section 5.1 we will assume N y ^> 1 and arbitrary D y , while in Section 5.2 
we will study the case with arbitrary N y and D y ^> 1. 

In general, the coherence length in the y direction (calculated at A9 X = 0, but 
trivially extendible to the case A9 X ^ 0), £ CJ/ is a function of D y , N y and z , as 
it can be concluded by inspection of the general expression for G in Eq. (67). 

As Ny ^> 1 and D y is arbitrary we will demonstrate that C, cy can be approxi- 
mated as 



n a 



(75) 



where $ is a universal function of the dimensionless parameters D y and A y . 
We will first study the asymptotic cases D y ^> 1 and D y 1, which are useful 
for physical understanding. Then we will generalize our results accounting for 
the influence of a finite divergence of the electron beam on the cross-spectral 
density. An analytical approximation for the function <3> will be proposed. This 
will be chosen to match in a very simple way both asymptotic expressions for 
D y ^> 1 and D y <^ 1: we will demonstrate that for any value of D y , dis- 
crepancies between the approximated and the actual (numerically calculated) 
expression are less than 10%, though there is no theoretical reason to assume, 
a priori, this relatively good accuracy. Also, at A y <^ 1 (and N y ^> 1) we will 
see that the VCZ theorem always hold. 

As N y is arbitrary and D y 3> 1 we will show, instead, that the coherence length 
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(calculated again at A9 X = 0, but trivially extendible to the case A9 X ^ 0) 
can be approximated as 



cy 



D y ,z 2 jN eS (N y ) 



(76) 



N e ff(Ny) being a universal function of the Fresnel number N y . \1/ is a universal 
function of the dimensionless parameters D y , N y and z . As usual in this 
paper, we will first study the asymptotic cases N y 1 and N y <^ 1. Then 
we will generalize our results accounting for any value of N y . As before, an 
analytical approximation for ^ will be proposed. This will be chosen to match 
in a very simple way both asymptotic expressions for N y ^> 1 and N y <^ 1: 
again, we will demonstrate that for any value of N y , discrepancies between the 
approximated and the actual (numerically calculated) expression are less than 
13% though there is no theoretical reason to assume, a priori, this relatively 
good accuracy. 

The case for a large Fresnel number with a finite electron beam divergence, 
or viceversa of a large beam divergence and a finite Fresnel number, is of 
practical importance. To give a numerical example, let us put A = 1 A, and 
consider a typical vertical emittance (for third generation sources in operation) 
e y ~ 3 ■ 1CT 11 m, that is e y ~ 2. On the one hand, if f3 y = 3, we have D y ~ 0.6 
and N y ~ 6. On the other hand, if f3 y = 0.3 we have, viceversa, D y ~ 6 and 
A> ~ 0.6. 



5.1 Very large Fresnel number N y ^> 1, arbitrary divergence parameter D, 



As a matter of fact, the only important assumptions used to derive results 
Section 4.1 were that N x ^> 1 and D x ^> 1 that are valid here as well. When 
these assumptions are granted, results can be factorized as a product of factors 
dependent separately on the x and on the y directions. Then the spectral 
degree of coherence in the horizontal direction will always be the same as in 
Section 4.1. Differences will arise here, of course, due to 9 yl ^ 9 y2 . 



5.1.1 Case with divergence parameter D y ^> 1. 

This case is the easiest to analyze, because one can follow step by step the 
derivation for the one dimensional model given in the previous Section 4.1. 
Calculations in the vertical direction y simply follow the derivation for the 
horizontal direction x. As a result, one can start from Eq. (67) and perform, 
separately for the x and the y directions, the same simplifications which hold 
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for the one-dimensional model. By comparison with Eq. (60), one obtains 
directly the result 



g(z Q , A9 X , A6 y ) = exp 



9 A D 

" O X 



A x + D x 



exp 



04 n 

y y , z 2 A9 2 y 



- Ay + Dy 



(77) 



Normalization of Eq. (77) according to Eq.(66) has been obtained simply 
setting G(z Q , 0, 0) = 1. In this case as well, two-dimensional Fourier Transform 
of \S[z Q , (f) x + (f)y}\ 2 calculated with respect the variables — 2A x z Q A9 x / (A x +D x ) 

and — 2A y z o A0y /(A y + D y ) gives, in analogy with Eq. (58), an unessential 
factor 



d(p x I d(p y S z , (f) x + 0, 



constant , 



(78) 



which has been included in the normalization. Again, similarly as before one 
can calculate the coherence area fl c {zo) defined in analogy with £ c (i ) as 



^c(^o) ^cx(Zo)£,cy(Zo) 



(79) 



Performing the integration yields: 



1/2 



(80) 



In the limit for a large value of z Q the coherence area exhibits quadratic de- 
pendence on z a , that is Q c — ► ■KZ 2 /(N x N y ) 1 / 2 while for z — > 1/2, that 
is at the end of the undulator, it converges to the constant value Q c — > 



7T 



1/2 



1/2 



l/(4Ay + 1/D X \ [l/(AN y ) + 1/Dy 
It should be noted as before that the asymptotic behavior for z Q 3> 1 of g 



g(z ,A6 x ,A6 y ) = exp 
and Q c 



-1A x z 2 AQ 2 x 



exp 



-2A y z 2 A6 2 y 



(81) 



a- nz ° 



(82) 



are direct application of van Cittert-Zernike theorem, as it must be since 
A XjV 1. In fact, Eq. (81) is simply linked with the two-dimensional Fourier 
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transform of F L (i x )F iy (l y ). 



5.1.2 Case with divergence parameter ^ < 1. 

With respect to the situation treated in Paragraph 5.1.1, this case requires a 
more careful analysis of the relations between different parameters. In fact, on 
the one hand D y <C 1 implies that the electron beam divergence drops out of 
the problem parameters, but on the other hand in this case the divergence of 
the radiation is described by the intrinsic divergence of undulator radiation, 
that is described by a more complicate mathematical function, compared with 
a Gaussian. In relation with this, it should be noted that simplifications in 
Section 4.1 pr 5.1.1 were based on the very specific properties of the Gaus- 
sian function representing the electron distributions in offset and deflection. 
Luckily, this is a realistic description in storage ring beam physics. 

Let us consider Eq.(74). In order to derive it we only used the assumptions 
N x ^> 1 and D x ^> 1, for 9 X = 9 y = 0. The result of operations on the right 
hand side of Eq. (74) depend on how N y scales with respect to z 2 . The cases for 
A y > 1 cannot be dealt with fully analytically. In the following we will analyze 
the asymptotic situation A y <^ 1 and then we will treat semi-analytically the 
generic situation for all values of A y . As we will see, as soon as A y < 1 we 
start to be in the applicability region of the VCZ theorem. 

It should be noted that the dependence in A9 X and A9 y in Eq. (74) are already 
separated. Therefore, what has been said in Section 4.1 regarding the behavior 
of coherence properties in the horizontal direction hold independently of the 
behavior of coherence properties in the vertical direction. 

(A) Far zone case: A y <C 1. — Since N y 1 it must be z% N y 1 in 
order to allow for A y <C 1. Since we are working in quasi-homogeneous source 
condition (N x ,y 3> 1) in the limiting situation A y <C 1 we should recover VCZ 
theorem: this is the far field case. Even in the presence of the extra parameter 
D y we can treat the generic case A y <C 1 independently on how A y compares 
with respect to D y . 

Eq. (74) can be simplified on the assumptions A y <C 1 and D y <C 1. In fact, 
the Gaussian exponential factor inside the integral in d(p y in Eq. (74) imposes 
a maximal value 0^ ~ A y + D y <C 1 . 

Simultaneously, from the oscillating factor, always inside the integral in d^> y , 
we have a condition for the maximal value of A9 y ~ < 1: in fact, if this 
condition is not fulfilled the integrand will be highly oscillatory. Alternatively, 
we can obtain a similar condition from the Gaussian exponential factor in A9 y 
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outside the integral, since A9 2 ~ \/{A y z 2 ) <C 1. 

As a result, the dependence of S on (0 + A9 y ) 2 can be dropped giving 



G{z ,Ad x ,A6 y ) =exp 
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The integral in d(f> y can be performed giving 
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G(z a , A9 X , AOy) — exp 
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Normalizing G in such a way that G(i , 0, 0) = 1 we obtain the following 
expression for the spectral degree of coherence g: 



g(z , A9 X , A9 y ) = exp 
x exp 
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Finally, combination of the second and the third exponential functions yields 
the result 



g(z , A9 X , A9 y ) = exp 



2A x z 2 A9 2 x D x 



(A x + D x ) 



exp 



-2A y z 2 A9 2 y 



(86) 



that is, again, a direct application of the van Cittert-Zernike theorem. For 
A6> = we have 
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(87) 
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(B) Case A y ^> 1. — This case encompasses situations with z Q 
well as situations with z Q ^$> 1. 



1 as 



Let us first consider z Q ^> 1. Looking at the integral in 4> y in Eq. (74) it is easy 
to recognize that its integrand is highly oscillatory in (p y when 2(p y z A9 y ^> 1, 
since A y j (A y +D y ) < 1. Therefore, the integrand will contribute to the integral 
significatively only up to values 2<f) y z A9 y < 1. On the other hand, the terms 
in S give non negligible contributions only for values of <j) y up to order unity. 
As a result, it must be 2z Q A9 y < 1. As z ~ 1 the width of g in A9 y is then 
of order unity. When z Q becomes larger than unity, the angular width A9 y 
will decrease and asymptotically, as z ^> 1, one will have a rapidly oscillating 
integrand for A9 y ~ l/z Q <C 1. However note that A9 y must be multiplied 
by z Q in order to obtain the correlation length which means that this remains 
constant and comparable with the diffraction length ^ cL w /uj as z Q increases. 

As a result, under the assumption z Q 1, terms in A9 y can be dropped in the 
functions S(-) of Eq. (74), and the functions S can be substituted with their 
limiting form sine. Moreover, in the limit for A y ^> 1 and D y C 1 we have 



G(z ,A9 x ,A9 y ) =exp 
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. exp 



A x + D x 

oo 

(2z A9 y ) 4] f d4> x sine 2 [(<g + 2 ) /a] , (88) 



where the simplification in the phase under the integral in d<f) y is possible 
for z A9 y (f) y D y / A y <^ 1 and the exponential function exp [—<j) 2 y /2 (A y + D y )\ 
under the integral in d<p y can be neglected because the sinc(-) function has 
characteristic length in 4> y of order unity. Eq. (88) is therefore valid in the 
limit for A y 1, D y 1 and z D 3> 1. The integral in d(j> y in Eq. (88) is 
simply the Fourier transform of the universal function 



htiy) = J d <t>x sine 2 [(4>l + 2 ) /4 



(89) 



done with respect to the variable 2z Q A9 y , conjugate to <f> y , that is 



G(z , A9 X , A9 y ) — exp 
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Fig. 15. The normalized radiant intensity Is, calculated from Eq. (89), as a function 
of the normalized vertical angle <f> y . 

It is not difficult to see that the width of this Fourier transform in A9 y z Q is 
much smaller than the characteristic width imposed by the Gaussian expo- 
nentials outside the integration sign, their ratio being of order D y <C 1. This 
means that the Gaussian in A9y outside the integral sign is almost constant 
with respect to the behavior of the Fourier transform, and can be neglected, 
to obtain 



G(z , A9 X , A9 y ) — exp 



2A x D x zjA9l 
A x + D x 



-1 00 


'<fi y exp 


J —00 





-i (2z A9 y ) 



The assumption z Q ^> 1 was vital for the derivation of Eq. (91). However, 
for all values of z Q such that A y 3> 1 (and, therefore, up to the exit of the 
undulator at z Q = 1/2), it is easy to see from Eq. (67) that the source is quasi- 
homogeneous, so that the "anti" VCZ theorem applies and the cross-spectral 
density of the field at the source plane forms a Fourier couple with the angular 
distribution of the radiant intensity, which is simply given by Is as defined in 
Eq. (89). 
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It is interesting to justify the fact that Is is in fact the angular distribution of 
the radiant intensity. To this purpose, it is sufficient to note that Is is simply, 
normalization factors aside, Eq. (176) of which represents the intensity 
from a beam with e x — > oo and e y — > 0. A representation of Is(4>y) is given in 
Fie. 15. 



As a result, one may therefore conclude that Eq. (91) is valid in general, for 
any value of z D such that A y ^> 1, in the asymptotic limit D y <C 1. 

It is important to note that Eq. (91) can be written as 



g(z a , A9 X , A9 y ) = exp 
where j(z A9 y ), given by 



2A x D x z 2 Q A9 2 x 



l{z A9 



y) > 



(92) 



j{Z O A0y 



2tt 2 



, cxp 



—i \ 2Z A6 y J 



(93) 



is a universal function normalized to unity. It is possible to calculate Eq. (93) 
analytically. To this purpose, it is sufficient to note that the Fourier Transform 



oo oo 



7i(C,v) = / # 



oo — oo 



exp 



sine 



2 ( + ^y 



can be evaluated with the help of the Bessel-Fourier formula as 



7! (A) = 2n J d(f) 0J O (<M) sine 2 (^j 



2vr 



Ti + A 2 Ci 



A 2 



— 2 sin 



A 2 



2Si 



(95) 



(94) 



where A 2 = £ 2 + 77 



2 A$ 



b\ + 2 , Si(-) is the sine integral function and Ci(-) 



is the cosine integral function. As a result one has 



j(Z A9y 
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(96) 



The function j(z A9 y ) is illustrated in Fig. 17. 
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Fig. 16. Values of exp[— 4> y /(2A y )]Is((j)y) as a function of 4> y . Solid line: A y = 10, 
the plot is still similar to Fig. 15. Dashed line: A y = 1.0. 



It should be noted that ^y(z A9 y ) is the spectral degree of coherence g calcu- 
lated for A9 X = in the limit for A y ^> 1 and D y <C 1, in agreement with Eq. 
(92). 



(C) Case A y 



In case (B) we have shown that terms in A9 y can 



be dropped in the functions S(-) of Eq. (74), for every value of z a , and the 
functions S can be substituted with their limiting form sine. In the case of 
A y ~ 1 and D y <C 1 we have 



G(z Q , A9 X , A6y) = exp 



x 



2A X D X %A6 2 X 
A x + D x 

t> y exp 



exp 



-2D y z 2 A9 2 y 
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Fig. 17. Behavior of | —y | as a function of z o A0 y . This universal plot illustrates the 
absolute value of the spectral degree of coherence \g\ calculated for A6 X = in the 



limit for A y S> 1 and D y <C 1 . 



The integral in d(f> y in Eq. (97) is the Fourier transform of 
Fs(4> y ) = Is(4>v)exp 



2 A,, 



(98) 



done with respect to the variable 2z Q A9 y , conjugate to 4> y . Similarly as before, 
it is not difficult to see that the width of this Fourier transform in A9 y z Q is 
much smaller than the characteristic width imposed by the Gaussian expo- 
nentials outside the integration sign, their ratio being of order D y <^ 1. This 
means that the Gaussian in A9 2 y outside the integral sign is almost constant 
with respect to the behavior of the Fourier transform, and can be neglected. 
Using the definition of F$ in Eq. (98) we have 



G(z , A9 X , A9 y ) — exp 



2A x D x z i Mj 
A x + D x 



,exp 



-i2z A9 y (f) y 



Fs(4> y ) 



(99) 
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Fig. 18. The behavior of the absolute value of the spectral degree of coherence g 
(at perfect resonance) as a function of z o A0 y . Here z Q = 1/2, N y = 10 and D y is 
negligible. The black circles represent actual numerical data. 

At this point, no simplification is possible and numerical analysis of the prob- 
lem should be undertaken in order to calculate G, followed by normalization 
according to 



G(z o ,0,0) = J 4 v F S (bv) , (100) 

— oo 

as already said, in order to find an expression for the complex degree of coher- 
ence g. This discussion underlines again how the correlation angle A6 y in the y 
direction behaves. It starts from a constant value equal to the diffraction angle 
\jcj (ujL w ) when z Q ~ 1, which corresponds to the maximal possible value of 
A y once N y is fixed. Then it decreases as z Q grows. Asymptotically in limit for 
z ^> 1 (but still such that z 2 Q < N y ), it behaves as ~ l/z Q , as it is clear from 
the fact that the function F$, which is Fourier Transformed in Eq. (99), does 
not depend on any parameter. However, it should be noted that A9 y must be 
multiplied by z Q in order to obtain the correlation length, which means that 
this remains constant and comparable with the diffraction length ^JcL w /u as 
z increases. Correlation length in the x direction is governed instead by the 
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Fig. 19. Normalized coherence length £ cy as a function of the normalized distance 
z when the electron beam divergence is negligible. Here N y = 10. The black circles 
represent actual numerical calculations. The asymptotic limit for A y <C 1 (VCZ 
theorem) is shown with a dotted line. Finally, the solid line is calculated with the 
approximated formula (101). 

Gaussian exponential function in A9 X outside the integral sign in Eq. (99), 
exactly as in the simplified model treated in Section 4.1. 

We performed some numerical calculation with the aim of giving the reader 
an exemplification. We set A6 X = and N y — 10. Assuming that we are in the 
asymptotic limit D < 1 we can rely on what has been said in this paragraph 
to calculate g and on Eq. (63) to calculate the correlation length £ cy in the 
vertical direction for any value of z . These numerical results must agree with 
the Van Cittert-Zernike limit for i s < 1 treated in paragraph 5.1.1 (A): in 
fact, Eq. (86) and Eq. (63) yield immediately a linear dependence of £ cy on 
z„ as A y < 1. Then, using paragraph 5.1.1 (B) we can extend the function 
£,cy{z ) to all values of z Q . In this way we obtain ^ C y{z ) for every value of z . 
In Fig. 16 we plot Eq. (98) for two particular values of A y : A y = 10 (solid 
line) and A y — 1 (dashed line). Note that the solid line Fig. 16 is still similar 
to Fig. 15, although the tails are changing and the width is already smaller. 
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The following step is to calculate the Fourier Transform of F$ according to 
Eq. (99) after normalization procedure according to Eq. (100). This gives g 
calculated for A9 X = 0. Fig. 18 illustrates | g | for z Q = 1/2 as a function of 
z o A0 y at A9 X = 0. Note that for this parameter choice we are still in the limit 
for A y >> 1, which explains why Fig. 18 is practically identical to the universal 
plot Fig. 17. 

One can calculate the coherence length i cy (z ) straightforwardly by means of 
Eq. (63). The curve obtained can be then compared, in the limit for A y <C 1, 
with the van Cittert-Zernike behavior illustrated in Paragraph 5.1.2 (A). Fig. 
19 shows our results. The black circles represent actual numerical calculations. 
The asymptotic limit for A y <C 1 (VCZ theorem) is shown with a dotted line. 
When N y >> 1 and D y <^ 1, the coherence length £ C y(z ) can be calculated 
with the approximated formula 



Scy — 



2 7F 

a 2 + — 
A 



n l/2 

(101) 



y-i 



Under the approximation of negligibly small electron beam divergence in the 
vertical direction, the normalized coherence length is thus a universal function 
of one dimensionless parameter, A y . On the one hand, Eq. (101) accounts for 
the asymptotic behavior as A y <C 1 (VCZ theorem). On the other hand, the 
value of £ cy in Eq. (101) approaches the constant value £ cy — > a for asymp- 
totically large values of A y . Beside accurately reproducing the asymptotes for 
small and large values of A y , Eq. (101) provides an accuracy of several per 
cent with respect to the result of numerical calculations (when A y is within 
the limits (0, oo)), in the whole range of the parameter A y . The solid line in 
Fig. 19 is calculated with the approximated formula (101), where calculated a 
numerically using Eq. (91) to calculate £ cy and obtaining a ~ 1.12. Therefore, 
according to Eq. (101), when N y ^$>l and D y <C 1, we have ^ cy = const = 1.12 
at the undulator exit (with accuracy D y <^ 1 and 1/N y 1) and, in dimen- 
sional units, £ cy = 1. \2\J L w c/uj. In this case the coherence length is a function 
of undulator length and wavelength due to the intrinsic divergence of the 
undulator radiation and a = 1.12 is a universal constant. 



5.1.3 Case with finite divergence parameter D y . 

We will first discuss, in Paragraph 5.1.3 (A), the limit for A y <^ 1: in this 
case, whatever the value of D y , we will recover the VCZ theorem. Note that, 
although we are discussing the limit iV" y ^> 1 there always be values of z Q large 
enough so that A y <^ 1. Further on, in Paragraph 5.1.3 (B), we will discuss 
the case A y 1: note that, since we are discussing the limit for N y ^> 1, this 
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will always be the case near the exit of the undulator at z — 1/2. We will be 
particularly interested to this situation; in fact, the study of the case A y ^> 1 
near the exit of the undulator will allow us to give an explicit representation 
of Eq. (75). 

(A) Far zone case A y <^ 1. — Eq. (74) is still valid and can be calcu- 
lated numerically, in principle, for any value of z . In the case A y <C 1 it can 
be further simplified. In fact, in this situation, the second exponential func- 
tion of the right hand side of Eq. (74) limits the possible values of A9 y to 
A0 y <C 1. Moreover z 3> 1 so that Eq. (97) is valid in this case. In particular, 
remembering the definition of F s (4> y ) in Eq. (98) we have 



G(z , A9 X , A9 y ) — exp 



2A x z 2 A9 2 x D x 



x 



(A x + D 

oo 

J d(j)y exp 



exp 



— - ( 

Ay + Dy 



2A y z 2 A9 2 y D y 



(Ay + Dy) 



(102) 



It is not difficult to see that the ratio between the characteristic width in 
A9 y of the exponential function outside the integral sign and the exponential 

function in inside the integral sign in Eq. (102) is of order \J A y j D y and it is 

always much smaller than unity unless D y <C 1: such a case has already been 
treated before in Paragraph 5.1.1 (A), and has been shown to obey the VCZ 
theorem 5 . For all other values of D y we have, automatically, D y ^> A y , so 
that we can neglect the integral in d$ y and we get back once more the Van 
Cittert-Zernike regime. To sum up, we obtain the following expression for g, 
which is valid for A v 1 with no restrictions on D v : 



g(z ,A9 x ,A9 y ) = exp 



2A x z 2 A9 2 x D x 
(A x + D x ) 



exp 



-2A y z 2 A9 2 y 



(103) 



Calculation of the coherence length from Eq. (103) at A9 a 
more the behavior 




1/2 



gives once 



(104) 



5 Alternatively one may note directly that <j) y can only range over values much 
smaller than unity. As a result, the dependence of Is on (f> y can be dropped, giving an 
extra normalization constant to be disposed of. Then, the integral in d(j) performed 
giving, as in Paragraph 5.1.1 (A), exp [—2A 2 z 2 A9 2 /(A y + D y )] to be combined with 
the exponential function in A6 y outside the integral sign, giving exp [— 2A y z 2 AQ 2 ]. 



60 



which is consistent with the partial result in Paragraph 5.1.1 (A). 



(B) Near zone, A y >> 1. — Equation (97) is still valid in this case and 
following the same line of reasoning as paragraph 5.1.2 (B) one gets 



g(z ,A9 x ,A6 y ) = exp 



2A x D x z 2 A6 



(A x + D x 



exp 



2A y D y z 2 A6 2 y 

(Ay+Dy) 



l{z A6 y ) , 



(105) 



valid for A y ^> 1 and arbitrary D y . Note that since we are working in the limit 
for N y 3> 1, for z Q = 1/2 we have A y 3> 1. We can see from Eq. (105) that, 
for finite values of D y , the cross-spectral density g, evaluated at the exit of 
the undulator for A8 X = 0, is given by the product of the exponential func- 
tion exp [-2D y z 2 A6l] with the funct ion illustrated in Fig. 18. This remark is 
intuitively sound. Since N y >> 1 in fact we have weakly quasi-homogeneous 
wavefronts near the exit of the undulator, and we can use the "anti" VCZ 
theorem to conclude that g must for a Fourier couple with the intensity distri- 
bution in the far zone. This will simply be, for any arbitrary D y , a convolution 

between a Gaussian distribution with rms width equal to \p^~ y and Is, which 
is the angular distribution of radiant intensity for a beam with e x — > oo and 
e y — > 0. Finally, the Fourier transform of such a convolution between two 
function is simply given by the product of the Fourier transforms of the two 
functions. 

(C) Approximate formula. — With in mind Eq. (80), Eq. (101) and Eq. 
(104) we make the working hypothesis that Eq. (75) has the form 



(71 71 
D eS (Dy) + Ay 



1/2 



(106) 



Within the assumption N y >> 1, we have seen that if D y 3> 1 Eq. (80) is valid 
with relative accuracy 1/D y . This means that, in this limit, D ef[ (Dy) = D y . 
Then we have seen that if D y <^ 1 Eq. (101) holds with accuracy 1/N y . This 
means that D c g(0) = 2.50. Moreover we have seen that in the far zone case, 
in the limit for A y 1, the VCZ theorem holds, in agreement with Eq. (104). 

We are now in position to calculate _D e fr = f(D y ) for any value of D y . Eval- 
uation of Eq. (105) at z Q = 1/2, followed by normalization according to 
G(5 o ,0,0) = 1 gives the function g. Further integration of \g\ 2 to calculate 
the correlation function allows to recover D eS = f(D y ) as plotted in Fig. 20. 

One may choose to calculate D c $ = f{D y ) numerically, but it is also possible 
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Fig. 20. Comparison between exact and interpolated D c g functions. Solid line: plot 
of the exact result. Circles: approximation according to Eq. (107). 

to use, with reasonable accuracy, the following analytical interpolation of D eS : 



D eS ~ D eS (0) + D y = 2.50 + t) y . 



(107) 



There is of some interest to compare the exact and interpolated D c g(D y ) 
functions. Fig. 21 shows the function 



A(Dy) = 



b eS (b y 



D cS {0) + Dy 



(108) 



There is seen to be good agreement between the interpolated and exact D c g 
functions for small and large value of D y . Noticeable discrepancies for D y close 
to unity are, anyway, less than 10%. 

Our conclusive result is, therefore, the following: when N y ^$>l and D y assumes 
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Fig. 21. Relative accuracy A(D y ) of the match between true and interpolated Z) e ff 
functions. 



arbitrary values we have: 



7T 



£ ~ I : + — 

CV V2.50 + D y Ay J 



\ 1/2 



(109) 



Also, it is important to remember that under conditions N y ^> 1 and arbitrary 
D y the spectral degree of coherence is given by Eq. (105) in the near zone, 
and by Eq. (103) in the far zone. 



The final step is to check that our main work hypothesis, i.e. that the coherence 
length has the form in Eq. (106) is correct. This can be done comparing Eq. 
(109) with numerical calculations for any given value of N y ^> 1 and a finite 
D y , which give a good agreement. 
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5.2 Very large divergence D y >> 1, arbitrary Fresnel number N y 



We now move on to treat the case with arbitrary beam transverse size com- 
pared with the diffraction size (i.e. arbitrary Fresnel number N y ) and large 
divergence compared with the diffraction angle (i.e. D y 3> 1). The partic- 
ular case for N y ^> 1 and D y ^> 1 overlaps with the previous Section 5.1 
and has already been treated in Section 5.1.1. The conclusion was that £ cy = 
(n/Dy + n / Ay) 1 / 2 . In all the other remaining cases Eq. (74) can be simplified 
as follows: 
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that is 



g(z , A9 X , A6 y ) = exp 



2A x z 2 A9lD x 
A x + D x 
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2A y z 2 A0 2 yDy 
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111) 
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As is shown in Appendix C, having defined 



oo oo 
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we have the important result 

f(z ,A9 y )=(3(A§ y ) 



nSmc 



K + [<Py 



Ad y f 



sine 



l 2 x + (^ + AOy) 2 



, (113) 



(114) 



for every choice of z D . (3 is defined in such a way to be normalized to unity. If 
we account for Eq. (114) we obtain the following expression for the spectral 
degree of coherence g: 
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Fig. 22. Plot of \/3\ as a function of A9 y . This universal plot illustrates the absolute 
value of the spectral degree of coherence \g\ calculated for A6 X = in the limit for 
N y < 1 and D y > 1. 



g(z , A9 X , A9 y ) = exp 



2A x z 2 A9 2 x D 3 



A x + D x 



exp 



2A y z 2 A9 2 y D y 

Ay + Dy 



(3(A9y) . 



(115) 



It should be noted that, as N y ^> 1, the width of the gaussian function in 
A9 y in Eq. (115) becomes much smaller than unity and the function (3 can be 
considered constant and drops out of the normalized expression for g. So, even 
if we did not analyze here the limit for N y ^> 1 (we did it in Paragraph 5.1.1), 
we see that the limit of Eq. (115) for N y ^> 1 restitutes the results found in 
Paragraph 5.1.1. Therefore we conclude that Eq. (115) is valid for any value 



of N, 



y 



As is shown in Appendix C, the function (3(A9 y ) can also be calculated as 



1 °r ( A9 " 

(5{A9 y ) = -J da a J l"^ 1 



7T - 2Si(« 2 



(116) 
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Fig. 23. Comparison between exact and approximated N e s(N y ) functions. Solid 
line: plot of the exact results. Circles: interpolation according to Eq. (121). 

where Si indicates the sine integral function. The representation of /3 in Eq. 
(116) is easier to deal with numerically, because it involves a one-dimensional 
integration only. Performing the integral, one can tabulate \/3\ to obtain the 
plot in Fig. 22. This is an universal plot. It should be noted that (3(A9 y ) is 
the spectral degree of coherence g calculated for A9 X = in the limit for 
N y <C 1 and D y ^> 1. Its generalization for arbitrary N y is given by Eq. (115). 
Using Eq. (115) and the tabulated values for the universal function (5 we can 
therefore calculate g numerically for any choice of N y and subsequently, we 
can calculate the coherence length £cy(z ). For instance, in the particular case 
N y — 10 and A9 X = we obtain the simple linear behavior 

Ccy(Zo) = 0.54S o (117) 



For a generic value of N y , one can introduce an effective function N e s(N y ) so 
that 



= \kkr (118) 
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On the one hand, the function N e s(N y ) can be computed numerically, similarly 
as we did for the particular case N y = 10. N eS (N y ) is represented by the solid 
line in Fig. 23. 

On the other hand, one may also use an interpolation for N eS (N y ). First, 
numerical calculations tell us that, in the particular case N y — > 0, we have 

^cy(Zo) = /^i • (H9) 

Second, as D y ^> 1 and N y ^> D y we have 




(120) 



The simpler interpolated formula which satisfies both asymptotes is therefore: 
N eS (N y ) ~N y + 0.35 , (121) 



The interpolation of N e s(N y ) is represented by black circles line in Fig. 23. 



There is of some interest to compare the exact and interpolated N e s(N y ) 
functions. Fig. 24 shows the function 



A(N y ) = 



N eS (N y 



N eS (0)+N y 



(122) 



There is seen to be good agreement between the interpolated and exact N e g 
functions for small and large value of N y . Noticeable discrepancies for N y close 
to unity are, anyway, less than 13%. 

Since in the case D y 3> 1 and N y ^> 1 we concluded that £ cy = (ir/D y + 
ix/ Ay) 1 ! 2 , we can formulate the hypothesis that, for D y 3> 1 and generic value 
of N y one has 



7X 



^cy(Zo) — ( - 



7X 



1/2 



0.35 + N v 



(123) 



The final step is to check that such hypothesis is correct. This can be done 
comparing Eq. (123) with numerical calculations for any given value of D y ^> 1 
and finite N y , which give a good agreement. 
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Fig. 24. Relative accuracy A(N y ) of the match between exact and interpolated N e g 
functions. 



5.3 Conditions for the source to be quasi-homogeneous 



Up to this moment we discussed, for simplicity, the case 9 X = 9 y = 0. We 
will now treat the generic case with arbitrary 9 X and 9 y . This discussion will 
reduce to the relation between the weakly quasi- homogeneous condition and 
our assumptions N x >> 1, D x >> 1 and either N y >> 1 or D y ^> 1. When 
we deal with weakly quasi- homogeneous wavefront, the results found for 9 X = 
9 y = have extended validity for generic values of 9 X and 9 y . As already said 
before though, we will find that the wavefronts are not always weakly quasi- 
homogeneous in the vertical y direction. In this case, previously found results 
are only valid in the particular case 9 y = 0. Let us consider the situation in 
more detail. We will start with Eq. (73), that may also be written as 



G(z , 9 X , 9 y , A9 X , A9 y ) = exp [i2 (9 x z A9 x + 9 y z A9 y 
r 91 + AA x z 2 A9 2 x D x + 4tA x 9 x z A9 x 
X6XP [ " 2(A X + D X ) 
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5. 3. 1 Very large divergence parameter D y ^> 1 

Because of the properties of the S function <fi y can only change of a quantity 
A<j) y ~ 1, otherwise the S functions will drop to zero. Then, since i) 9 > 1 we 
have from Eq. (124) 
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J d$ x S*[z J 2 x + ($ y -A9 y ) 2 ]s[z J 2 x + ^ y + A9 y ) 2 } .(125) 



that is obviously weakly quasi-homogeneous. 

5.3.2 Very large Fresnel number N y >> 1. 

Case A y ^> 1. — In the case i 9 > 1 we can follow the same line of 
reasoning in Section 5.3.1 simply replacing the roles of D y with A y , obtaining 



G(z , 9 X , B y , AB X , AB y ) = exp [i2 (B x z AB x + 9 y z AB y ) 



x exp 



01 + AA x z 2 A9 x D x + UA x B x z a AB x 



exp 



2{A X + D X ) 

2AyZ 2 A9 2 Dy + 2iAy9yZ A9y 
Ay + Dy 



69 



x exp 
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2i<p y A y z A9 y 



x 



(Ay+Dy) J 

dfcS* z Q , 4> 2 X + (fa - A^) 2 S z Q , 2 + (fa + A6 y ) 2 ] . (126) 



that is obviously weakly quasi-homogeneous. 

(B) Case A y <^ 1. — In this situation it must be z > 1. Therefore we 
can substitute all S functions with sine functions. The case D y 1 has been 
already treated. Let us, therefore, first assume D y ~ 1. Using the fact that 
A9 2 < 1/z 2 < 1, Eq. (124) gives directly 



G(z , 9 X , 6 y , A9 X , A9 y ) = exp [i2 (9 x z A9 x + 9 y z A9 y ) 



x exp 



# 2 + AA x z 2 A9 2 x D x + ^AA^Afl, 



2 (A, + D x ) 



x 



exp 



exp 
(4>y + ^) 2 



-2A y z 2 o A0 2 y 



2A, 



(127) 



that is obviously weakly quasi-homogeneous. 

Now let us consider the case D y <C I. Eq. (124) can be simplified on the 
assumptions A y <C 1 and D y <C 1. In fact, the Gaussian exponential factor 
inside the integral in d(f> y in Eq. (124) imposes a maximal value (9 y + <p y ) 2 ~ 
A y + D y 1. 

Simultaneously, from the Gaussian exponential factor in A9 y outside the in- 
tegral, we have a condition for the maximal value of A9 2 ~ 1/ ' z 2 <C 1 since 
A9 2 ~ l/(i,4 2 ) « 1. 

As a result, the dependence of 5 on (0 + A9 y ) 2 can be substituted by a depen- 
dence on (9^ and, as has already been said, the S functions can be substituted 
with sine functions, giving 



G(z , 9 X , 9 y , A9 X , A9 y ) = exp [i2 (9 x z A9 x + 9 y z A9 



x exp 
x exp 



9j + 4A x zjA6 2 x D x + 4iA x 8 x z A6 x 

2{A X + D X ) 
2A y z 2 A8 2 D y + 2tA y 8yZ A8 y '] 

Ay + Dy 
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X 



j d(j)y 



cxp 



{9 y + 4>yf + 2<f>y (2iA y z A6 l 
2(4 + Dy) 



.(128) 



The integral in d(p y can be performed giving 



G(z , 9 X , 9y, A9 X , A9y) — exp %2 (9 X z Q A9 x + 9 y z A9 y 
9\ + AA x z 2 A9 2 x D x + AtA x 9 x z A9 : 



x exp 



x exp 



2{A X + D x ) 
2A y z 2 A9 2 y D y + 2iA y 9yZ A9 y 



x exp 



2A 2 yZ 2 A9 2 y 
Ay + Dy 



cxp 



Ay + Dy 
2iAy9yZ A9y 



Ay + Dy 



Is{9 y ) . (129) 



Normalizing G in such a way that G(z D , 9 X , 9 y , 0, 0) = 1 we obtain the following 
expression for the spectral degree of coherence g: 



g( Zo , A9 X , A9 y ) = exp [i2 (9 x z A9 x + 9 y z A9 y ) 



x exp 



9 2 x + AA x z 2 A9 2 x D x + AiA x 9 x z Q A9 x 
2(A X + D X ) 



cxp 



-2A y z 2 A9 2 y 



, (130) 



which generalizes Eq. (85) for any value of 9 X and 9 y , and shows weak quasi- 
homogeneity of the wavefronts. 

(C) Case A y ~ 1. — Since N y ^> 1 and z Q ^> 1, the exponential function 
in A9 y outside the integral sign of Eq. (73) impose A9 y <C 1 so that 



G(z , 9 X , 9 y , A9 X , A9 y ) = exp [i2 [9 x z A9 x + 9 y z A9 y 
9\ + AA x z 2 A9 2 x D x + 4iA x 9 x z A9 : 



x exp 



cxp 



2(A X + D x ) 

2AyZ 2 A9 2 y Dy + 2lAy9yZ A9y 
Ay + Dy 



x 



j d(f)y 



cxp 



(9y + (j>y) 2 + A4yAyZ A9y 

2(Ay + b y ) 



Istfy) ■ (131) 



The intensity distribution can be found from Eq. (131) setting A9 X = A9 y = 
thus obtaining 
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I(z , 9 X , Qy) = exp 



01 



[ 2{A X + D X )\ 



J d<j> y exp 



(Oy + kf 

2(Ay + Dy) 



(132) 



It is evident by inspection that we cannot factorize Eq. (131) to obtain \G\ = 
I(z , 9 X , 9y)w(z , A9 X , A9 y ) (where have put \g\ = w(z a , A9 X , A9 y )). As a result 
we conclude that, in this case, the wavefront is not quasi-homogeneous, not 
even in the weak sense. 



5. 3. 3 Discussion. 

Let us discuss the results obtained in this analysis of quasi-homogeneity. 

In Section 4.1 we have seen that weakly quasi- homogeneous wavefronts which 
are not quasi- homogeneous in the usual sense are present in the far field, when 
the VCZ theorem holds. From our analysis in the one-dimensional framework, 
the notion of far zone arises in the x direction, when the apparent angular 
dimension A x of the source is much smaller than the divergence of the radiation 
beam that, in this case, can be identified with the electron beam divergence 
4 > 1. Both N x >■ 1 and ^ > 1. If i x < £) x the wavefront is quasi- 
homogenous but only in the weak sense, one is in the far zone and the VCZ 
theorem applies. If A x ^> D x one is in the near field zone and the wavefront 
is quasi-homogeneous in the usual sense. Transition from the near to the far 
zone always involves, in this case, weakly quasi- homogeneous wavefront. 

In the two-dimensional framework studied in the present Section, as for the 
one-dimensional model, both N x and D x ^> 1. The x and the y coordinates 
appear factorized and the far zone applies now separately to both x and y 
directions, meaning that A X)V of the source is much smaller than the divergence 
of the radiation beam. In other words, we are in the far zone as soon as 
the (square of the) beam size at a given z Q begins to be much larger than 
the (square of the) initial size of radiation (i.e. much larger than the Fresnel 
numbers N Xjy ). 

In the case D y >> 1, independently on the value of N y , we will always have 
weakly quasi- homogeneous (but not always quasi-homogeneous in the usual 
sense!) wavefronts at any distance z Q . This is exactly the situation discussed 
for the x direction, where transitions from the near (A y ^> D y ) to the far field 
{A y Dy) involve only weakly quasi- homogeneous wavefronts. Moreover, in 
the far field zone, the VCZ theorem is valid. 

In the case with D y < 1, N y ^> 1 we will have quasi- homogeneous wavefronts in 
the usual sense at z Q ~ 1 (near zone) and weak quasi-homogeneous wavefronts 
in the far zone, at z a ^> 1. 
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Finally, if N y 3> 1, D y < 1 and A y ~ 1, i.e. for distances z ~ yiVy, we have 
seen that the wavefront is non quasi-homogenous, not even in the weak sense. 

These results can be seen in terms of convolution between the Gaussian dis- 
tribution of intensity associated to the electron beam emittance and the dis- 
tribution of intensity due to intrinsic properties of undulator radiation. The 
only case, among those treated up to now, when such a convolution does not 
simplify into the product of separate factors is when N y 1, D y not much 
larger than unity and A y ~ 1: in this case we do not have quasi-homogeneous 
wavefronts, not even in the weak sense. 



6 Radiation from some non-homogeneous undulator sources 

In Section 4.1 we treated a simplified situation with N x ^> 1, D x ^> 1 and 
y \ = y2 . Moreover, for simplicity of calculations we assumed N y <^ 1 and 
D y <C 1. In Section 5, instead, we treated the case of electron beams with 
N x >> 1, D x >> 1 and either N y 3> 1 or D y 3> 1 (or both). It is important to 
note that, in the x direction, the results obtained in Section 4.1 are the same 
as the one in Section 5. In fact as N x ^> 1 and D x ^> 1, the cross-spectral 
density factorizes in the product of two contributions depending separately on 
the x and y coordinates, and under N x ^> 1 and D x ^> 1 the derivation of the 
x-dependent factor is always the same. 

We have seen that in some of the cases discussed in Section 5.3, the assump- 
tions N x ^> 1, D x ^> 1 and either N y ^> 1 or D y ^> 1 (or both) were enough to 
guarantee that the wavefront is weakly quasi-homogeneous in the sense speci- 
fied by Eq. (23). In this Section we will extend our analytical investigations to 
some cases outside the range of parameters treated before, where the weakly 
quasi-homogeneous assumption is not fulfilled in the far zone. In particular, 
we will demonstrate that, under conditions N x ^> 1, D x ^> 1 and both N y <C 1 
and < 1, wavefronts are not weakly quasi-homogenous in the vertical y 
direction (although they are in the horizontal x direction). 

First, in Section 6.1 we will analyze the case N x ~ 1, D x <^ 1. Assuming 
a vertical emittance of the electron beam much smaller than the horizontal 
emittance e y e x , we have, automatically ^ C 1 and ^ C 1. This cor- 
responds to a practically important situation. For instance, consider a VUV 
beamline at a third generation light source with A = 30 nm, e x = 3 • 10~ 9 m, 
and t y = 0.03 • 10~ 9 m, i.e. e x = 0.6 and e y = 6 • 10~ 3 . If j3 x — 3 we would 
have N x — 2 and D x = 0.2. Second, in Section 6.2 we will study the situation 
N y < 1 and D y 1 with N x 1 and 4 » 1, that will give us back also 
the limiting case for N y <^ 1 and D y <^ 1 with N x 3> 1 and D x 3> 1 (already 
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discussed in Section 4.1). The situation with finite vertical Fresnel number and 
negligible vertical divergence (compared with the diffraction angle) is a very 
practical one: for instance, given a third generation light source with A = 1A 
and e y = 1CT 11 m, i.e. e y = 0.6, a value (5 y — 6 corresponds to D y — 0.1 and 
N y = 3.6. 

Although in these two cases, the weakly quasi-homogeneous assumption is not 
fulfilled we will see that the choice z a 3> 1 will allow us to treat these situations 
in analogy with respect to some weakly quasi-homogeneous case we already 
dealt with. 



6.1 Case N x ~ I, D x < 1. 



Assuming a vertical emittance of the ring much smaller than the horizontal 
emittance e y <C e x , we have, automatically, N y <C 1 and D y <^ 1. We start 
with Eq. (67), that can be specialized to an equation dependent on the x 
coordinates only and, in the case for z Q ^> 1, can be written as 



G(z ,e,A9) 



exp (i20Af?£ o ) 

oo 

x J d(j) 



exp 



exp 



4NA6 2 D + 4i{N/z )6A6 
2{N/zl + D) . 

Ai4>{N/z o )A0 



—\ 2 

e 



2(N/% + D) 



xsmc 



((f) - Ad) 2 / A sine ((f) + AO) 2 /A , (133) 



where we systematically omitted x subscripts. Since A + D 1, the expo- 
nential factor in (0 + 0) 2 inside the integral sign in Eq. (133) behaves like a 
5-Dirac function with respect to the sinc(-) functions inside the same integral: 
as a result we can substitute <p with 9 in the sinc(-) functions, which drop out 
of the integral sign. Then, the integral in d(f> can be calculated analytically so 
that we have: 



G(z , 9, AO) = exp (i20A0z o ) exp 
(0- A9) 2 /A 



xsmc 



sine 



2NA9 2 D ' 
'n/z 2 + D_ 
9 + AO) 2 /A 



x exp 



2(N 2 /z 2 )A9 2 



■ (134) 



N/z 2 + D 

Combination of the second and the third exponential function yields 
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G(z , 9, AO) = exp (i29A9z Q ) exp [-2NA9 2 



xsmc 



(6 - AO) 2 /A sine (9 + A9) 2 /A . (135) 



Finally, in order to obtain the degree of spectral coherence g we should nor- 
malize G according to Eq. (66) 6 thus obtaining 



g(z , 9, A9) = exp (i29A9z Q ) exp [-2NA9 2 ] . (136) 

According to Eq. (136) the spectral degree of coherence g is such that \g\ is 
only function of A9. Moreover, the dependence of g on the phase 29A9z Q is a 
feature for radiation from Schell's model sources in the far field, that we have 
already encountered many times in the study of weakly quasi-homogenous 
cases. As a result we can conclude that the radiation of the undulator source 
at N x ~ 1 and D x <^ 1 represents the far field radiation of a Schell's model 
source. Finally, it should be noted that in a two-pinhole experiment, for any 
vertical position of the pinholes, the fringe visibility, i.e. the modulus of the 
spectral degree of coherence depends only on the separation along the hori- 
zontal x direction. As a result, while in Section 4.1 we put 9 y \ = 9 y 2, thus 
selecting from the very beginning a horizontal plane, the present case can be 
fully described by a one-dimensional model, independently on the choice of 
transverse coordinates of the pinholes. 



6.2 Case N y <\ and D y < 1 with N x > I, D x > 1 



In this situation we go back to a two-dimensional model. When 2„ > 1 we 
have A y 1, which is the limiting case treated in Paragraph 5.1.2 (A): the 
difference is that, now N y ~ 1. Let us start with Eq. (73) written as 



G(z , 9 X , 9 y , A9 X , A9 y ) = exp [i2 (9 x z A9 x + 9 y z A9 y ) 



x exp 



9 2 x + AA x z 2 Q A9 2 x D x + AiAj x z A9 x 



exp 



,cxp 



2(A X + D x ) 

2A y Z 2 A9 2 y by + 2iAy9yZ A9y 
Ay + Dy 

(9y + (/)y) 2 + At(j)yA y Z A9y 



2(4/ + Dy) 



6 Note that, in this case, we are dealing with non quasi-homogeneous wavefronts 
and, as has already been said, normalizing according to G(z o ,0, 0) = 1 is not the 
same of normalizing according to Eq. (66). 
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z , <t>i + (0, - Ae y y s z , <\>i + (0, + Ae y y , (137) 



A y <C 1 and D y <C 1 impose a maximal value of (9 y + 4> y ) 2 ~ A y + D y <C 1. 
Moreover the S 1 functions can be substituted with sine functions since we are 
working in the limit for z Q ^> 1. Then, from Eq. (137) we obtain 
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, (138) 



As done before, the integral in d<f> y can be performed giving 



G(z , 9 X , 9 y , A9 X , A9 y ) = exp [i2 (9 x z A9 x + 9 y z a A9 y 
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that is 



G(z , 9 X , 9 y , A9 X , A9 y ) = exp [i2 (9 x z A9 x + 9 y z Q A9 y 
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Finally, normalization according to Eq. (66) yields 



g(z , 9 X , 9 y , A9 X , A9 y ) = exp [i2 (9 x z A9 x + 9 y z A9 
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Obviously |p| is a function of both 9 y and A^, so that in this case we have 
neither weak quasi-homogeneity, neither wavefronts which can be described 
by Schell's model. However it should be noted that the integrals in Eq. (141) 
do not contain any parametric dependence. 

It is interesting to have some final comment on Eq. (141). In the limit for 
N y <^ 1, we recover the case discussed in Section 4.1, the only difference 
being that we did not set 9 y \ = 9 y 2 = 0: in Section 4.1 we chose to deal 
with a one-dimensional model putting ourselves on the horizontal plane. Now, 
Eq. (141) allows to study the full two-dimensional situation in the limit for 
N y 1. In this case we see that the exponential function exp [— 2A y z 2 A9 2 y ] 

can be neglected because we have a maximum value of A9 y ~ 1. As a result 
we obtain: 



g(z a , 9 X , 9 y , A9 X , A9 y ) = exp i29 x z Q A9, 



x exp 
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Fig. 25. Three-dimensional representation of % as a function of 9 y and A8 y . 

It should be noted that Eq. (143) does not depend on parameters and is, in fact, 
a universal function. Besides a geometrical factor exp [i29 y z A9 y ] , the function 
X represents the spectral degree of coherence in the vertical direction, once the 
horizontal coordinates are fixed. The fact that it is a universal function means 
that even in the case of zero vertical emittance we never have full coherence 
in the vertical direction. On the one hand, this phenomenon can be seen to be 
an influence of the presence of horizontal emittance on the vertical coherence 
properties of the photon beam, as the integral in d(f) x in \ comes from an 
integration over the horizontal electron beam distribution. On the other hand, 
being x a universal function, the influence of the horizontal emittance on the 
vertical coherence does not depend, in the limit for N x 3> 1 and Dj. > 1, on 
the actual values of N x and D x . It is straightforward to see that x is symmetric 
with respect to A9 y and with respect to the exchange of A6 y with 9 y . When 
y = 0, i.e. 9 y i = —0y2, we obviously obtain x(0, A9 y ) = 1 that corresponds to 
complete coherence. In Fig. 25 we plot the three-dimensional representation 
of x{dy,A9 y ). In order to get a feeling for the behavior of x we a ^ so plot, in 
Fig. 26 and Fig. 27, two cuts of Fig. 25 illustrating, respectively, the behavior 
of x for a fixed 9 y = 0.5 (or fixed A9 y = 0.5) and at 9 y = A9 y . 

As it is evident from Fig. 25, x exhibits, for any fixed value of A9 y , many 
different zeros in 9 y . In Fig. 28 we illustrate some of these zeros as a function 
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X(0.5, x) 



0,8- 



0,4- 



-0,4- 



-0,8- 




Fig. 26. Plot x(0.5, x), illustrating the cut of Fig. 25 at 6 y = 0.5 (or fixed A9 y = 0.5). 

of A9 y , 9y t z(A9 y ). The interest of this plot is that, once a certain distance 
z A9 y between two pinholes is fixed, it illustrates at what position of the 
pinhole system, 9 y t z, the spectral degree of coherence drops from unity to zero 
for the first time. 

It is interesting to compare Fig. 28, with the directivity diagram of the radiant 
intensity Is(9 y ). This comparison is shown in Fig. 29. In the limit for N y <^ 1 
and D y <C 1 , one may increase the degree of coherence of the beam by spatially 
filtering the radiation in the far field. If a vertical slit is used with aperture d 
much larger than the horizontal coherence length, i.e. d ^> l; cx , one would have 
poor coherence. Decreasing the aperture of the slit will increase the coherence 
of the X-ray beam up to some value d smaller than £ c . Within our assumption 
N x 3> 1 one has the far field approximation £ cx = (ir / 'Nx) 1 / 2 ^. When d 
becomes smaller and smaller with respect to (n / 'N x ) l l 2 z the spectral degree 
of coherence g can be identified with the universal function x, as once can see 
by inspecting Eq. (142). As a result, as d becomes smaller one loses photons, 
but the X-ray beam transverse coherence ceases to improve because, as is seen 
in Fig. 29, the transverse degree of coherence g = x drops to zero along the 
vertical radiation pattern of the filtered X-ray beam: for instance, from Fig. 
29 one can see that x drops to zero for the first time at A# ~ 1 9 y ~ 2, 
where the X-ray flux is still intense. This behavior of the degree of coherence 
should be taken into account at the stage of planning experiments. To give 
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X( x,x) 




Fig. 27. Plot x(x,x), illustrating the cut of Fig. 25 at y = A9 y . 

an example, after spatial filtering, one may conduct a two-pinhole experiment 
(like the one illustrated in Fig. 7) and find, surprisingly, that for some vertical 
position 9 y of the pinholes (at fixed A9 y ) well within the radiation pattern 
diagram he will have no fringes, but for some other vertical position he can 
find perfect visibility. So, without the knowledge of the function x a user would 
not even have the possibility to predict the outcomes of a simple two-pinhole 
experiment. 

From the definitions of x, P, 7 and I s it can be seen that all universal functions 
introduced in this work are partial cases of the more generic 
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Fig. 28. Plot of 9 V: z as a function of A0 y . are some zeros of x, i.e. some of the 
values of 0^ such that x(A0 y , = 0. 

J 5 (0 S )=M(0 S ,O) 

oo 

— OO 
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7(^) = 2^2 / ^expH(2x)e]M(e,0) .(145) 
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The knowledge of M is all one needs to calculate coherence properties out of 
many experimental setups, in very practical situations. It is therefore worth 
to tabulate M. We present a 3D plot of M, obtained from such tabulation, in 
Fig. 30. 

Finally, it is interesting to sum up and compare results for the far field region 
obtained in this Section (non-homogeneous undulator source) with results ob- 
tained in Section 5. Many users performing coherent experiments with X-ray 
beams are interested in the beam coherence properties in the far field. We have 
seen that in the most general situation for third generation light sources one 
is interested in the case N x ^> 1 and D x ^> 1, which guarantees factorization 
of results in the x and y direction, with arbitrary N y and D y . In this case we 
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Fig. 29. Comparison between some zeros of x, A9 z (9 y ) (black circles), and the di- 
rectivity diagram of undulator radiation in the vertical direction at very large hori- 
zontal electron beam divergence D x S> 1 and negligible vertical divergence D y <C 1 
(solid line). 

will not have, in general, weakly quasi-homogeneous radiation in the vertical 
direction. The spectral degree of coherence can be found by simplifying Eq. 
(73) in the mathematical limit z — > oo: 



g(z , 9 X , 6 y , A9 X , A9 y ) = exp [i2 (9 x z A9 x + 9 y z A9 y ) 
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Is(4>y)\ -(146) 



We can see that for any value of N y and D y , in the far field limit we obtain 
a contribution to the cross-spectral density for the x and for the y direction. 
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A0 y axis 



Fig. 30. Three-dimensional representation of M as a function of 6 y and A6 y . 

The contribution for the y direction can be expressed in terms of the product 
of an exponential function and convolutions between the (Gaussian) electron 
beam divergence and universal functions. 

On the one hand, as D y <C 1 we have 



g(z , 6 X , 6 y , A6 X , A6 y ) = exp %2 (6 x z A6 x + 6 y z A6 : 



x exp 



-2N X A9 



exp 



-2N y A0 2 y 



X(0 V ,M V ) .(147) 



On the other hand, as D y ^> 1 we have weakly quasi- homogeneous wavefronts 
and 



g(z Q , 6 X , 9 y , A6 X , A9 y ) = exp i2 (0 x z o A9 x + 9 y z A6, 



x exp 



-2N X A9 2 X 



exp 



2N y A§l P(A9 y ) (148) 



where (3{A6 y ) is given in Eq. (116). 

Moreover, as N y 3> 1, and for arbitrary D y we have: 
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g(z , 9 X , 6 y , A9 X , A9 y ) = exp i2 (9 x z A9 x + 9 y z A9 y ^j 



x exp 



-2N T A9t 



exp 



-2N y A9 2 y 



. (149) 



It should be noted that Eq. (149) is simply a consequence of the application 
of the VCZ theorem in both horizontal and vertical directions. 



7 Application: Coherent X-ray beam expander scheme 



In this Section we show how transverse coherence properties of an X-ray beam 
can be manipulated to obtain a larger coherent spot-size on a sample. 

The idea of increasing the horizontal width of the coherence spot is based 
on the use of a downstream slit for selection of the transversely coherent 
fraction of undulator radiation. Imagine a slit very close to the exit of the 
undulator with an aperture d comparable with the coherent length of the 

radiation at the exit of the undulator £ cx = \Jtt/D x , as illustrated in Fig. 31. 
The new radiation source after the slit is now coherent and characterized by 

a horizontal dimension of the light spot equal to \J ir/D x . In the far field one 
can take advantage of the reciprocal width relations of Fourier transform pairs 
or, equivalently, the expression for the Fraunhoffer diffraction pattern from a 
slit, i.e. a sine function, to calculate the magnitude of the coherence spot. 
There is, of course, some arbitrary convention to agree upon when it comes to 
the definition of the width of the sine function but, numerical factors aside, 

this reasoning shows qualitatively that the coherence spot is of order \[t) x z 
which is e x times larger than the spot size dimension in the case of free- 
space propagation, of order z j \[ff x . The radiation beyond the slit must be 
then spectrally filtered by a monochromator (not shown in Fig. 31) to further 
narrow the spectral bandwidth. Here we assume that the radiation frequency 
uj is equal to the fundamental frequency u Q . The radiation beyond the slit is 
transversely coherent when the aperture d is equal (at most) to the coherence 
length £ cx . 

Let us present a numerical example illustrating the improvement of the hori- 
zontal coherence length obtained by slit application. Let us consider the case 
when the electron horizontal emittance is large e x ^> 1 and the vertical emit- 
tance is small e y <C 1, with N x ^> D x . Since N x = i x /3 x and D x = e x //3 x , this 
is the case, for instance when /3 X ~ 10. Then, for e x = 100 we have N x = 1000 
and D x = 10. This particular numerical example has been considered in Sec- 
tion 4.1 to illustrate, in free space, the behavior of the coherence length as 
a function of the position along the beamline, given in Eq. (64). Suppose we 

install the slit at z Q = 2. From Eq. (64) we have ^ cx ^ o=2 ~ \Jn/ D x ~ 0.6. 
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for x - direction 



X c>Yc 





vertical 
slit 



transversely 
coherent 
X-ray beam 

Fig. 31. Undulator radiation with a coherent X-ray beam expander. 

At z a = 12 we have a situation close to the asymptotic behavior where the 
field diffracted by the slit can be treated in the Fraunhofer approximation and 

£ cx \z =i2 ~ \flXcZ ~ 30. Comparison of the asymptotic behaviors after spatial 
filtering with respect to the free-space propagation case is given in Fig. 32. 
Note that only the asymptotic behaviors near the slit and at large distance 
z Q 1 are plotted for the spatially filtered radiation. What is important is, in 
fact, the comparison between the coherent distance at large values of z Q with 
and without spatial filtering in the near zone. 

This is an example in which evolution of transverse coherence through the 
beam line plays an important role. In fact, the ability of spatially filter ra- 
diation by a slit requires the knowledge of the transverse coherence length 
variation along the beamline. 

Let us calculate the number of coherent photons observed beyond the aperture. 
In the region of parameters where N x ^> 1 and D x >> 1, the number of 



transversely coherent photons into the slit aperture d = £ c 

as 



(N ph ) 



dK 



ph ) coh 



ph 



dx 



can be calculated 



(150) 



In the near-zone limit the slit is positioned at a position down the beamline 
z s ~ j3 x so that, from Eq. (64) and Eq. (50) we have: 

dN^ = JU (151) 
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Fig. 32. Behavior of the horizontal coherence length as a function of distance for 
free space (dashed line) and comparison with asymptotic behaviors near the slit and 
at large distances (solid lines). 



and 



7T- 



2na 7 



(152) 



Therefore we can write 



dN ph XN ph 



£ cx = y ft-=^ = . (153) 



At the opposite extreme, with the distance z s much larger than (3 X we find 
that, always from Eq. (64) and Eq. (50) 

dx yjlvzlal 
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(155) 



Therefore we can write 



dN ph 



XN ph 



(156) 



dx 



v 7 ^ 




Thus, the number of transversely coherent photons into the slit aperture d = 
£ cx will be independent of the distance z s . This means that operation of spatial 
filtering in the near field region, as proposed by us, will not diminish the 
number of coherent photons with respect to the usual practice in which spatial 
filtering to obtain coherent radiation is performed in the far field: its effect will 
only be that of increasing the horizontal dimension of the coherent spot size. 
This possibility to create transversely coherent radiation with large divergence 
in the horizontal direction is important for many experiments. The distance 
of the slit from the exit of the undulator sets a limit to achievable linear 
dimension of coherence area. If we build a coherent X-ray beam line, we want 
to have large linear dimension of coherence area at the specimen position. 
Therefore, in order to have a largest linear dimension of coherence area at the 
specimen position ^ C x{ z o) — o~x /Z o we must have a slit aperture of at most of 
the size d ~ \/{2na x i) installed in the near zone at z s ~ f3 x . In order not 
to loose coherent photons instead, the slit aperture must be at least of the 
size d ~ A/(27r<7 x /). The right compromise is thus a slit aperture of the size 
d ~ \/(2na x >) installed in the near zone at z s ~ f3 x . 



8 Conclusions 

Before this work, no satisfactory theory describing spatial coherence from 
undulator radiation sources has been built. In this paper we developed such 
a theory of transverse coherence dealing with X-ray beams, with particular 
attention to third generation light sources. 

First we studied Synchrotron Radiation as a random statistical process using 
the language of Statistical Optics. Statistical Optics developed around Gaus- 
sian, stationary processes characterized by quasi-homogeneous sources; under 
these assumptions, the characterization of statistical properties of the process 
are greatly simplified and the van Cittert-Zernike theorem (or its generalized 
version) can be used in order to describe the X-ray beam partial coherence 
properties in the far field region. However, for Synchrotron Radiation, there 
is no a priori reason to hold these assumptions satisfied. 
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We showed that Synchrotron Radiation is a Gaussian random process. As a 
result, statistical properties of Synchrotron Radiation are described satisfac- 
tory by second-order field correlation functions. We used a frequency domain 
analysis to describe them from a mathematical viewpoint. This choice is very 
natural. In fact, up-to-date detectors are limited to about 100 ps time reso- 
lution: therefore, in real-life experiments with third generation light sources, 
detectors are by no means able to resolve a single X-ray pulse in time domain 
and work, instead, by counting the number of photons at a certain frequency 
over an integration time longer than the radiation pulse. As a consequence 
of the frequency domain analysis we could study the spatial correlation for a 
given frequency content using the cross-spectral density of the system which, 
independently of the spectral correlation function, can be used to extract use- 
ful information even if the process is not stationary. 

We gave an expression for the process cross-spectral density dependent on 
six dimensionless parameters. Subsequently we tuned parameters at perfect 
resonance, thus obtaining a simplified expression. 

First we studied the limit of applicability of the quasi-homogeneous model 
from an analytical viewpoint, within the framework of simplifying assump- 
tions, namely in the limit of small electron beam divergence and Fresnel num- 
ber in the vertical direction, of large electron beam divergence and Fresnel 
number in the horizontal direction, and performing calculations for the cross- 
spectral density on the horizontal plane only. This simplified study allowed us 
to introduce the concept of weakly quasi-homogeneous radiation and virtual 
quasi- homogeneous source while discussing the applicability region of the VCZ 
theorem. 

Second, we studied the effect of the vertical emittance on the cross-spectral 
density. This study led us to analyze both cases in which the sources are 
weakly quasi-homogeneous and cases when they are not quasi-homogeneous 
at all. In the limit for large horizontal beam divergence and Fresnel number, 
which is always satisfied for third generation light sources, we found that the 
spectral degree of coherence factorizes in the product of factors depending 
separately on the horizontal and on the vertical coordinates. In the far field 
limit the vertical part of the spectral degree of coherence can be expressed in 
terms of the product of an exponential function (which, alone, would simply 
satisfy the VCZ theorem) and convolutions between the electron beam diver- 
gence in the vertical direction and a universal functions, that we introduced 
in our work. The universality of such a function implied that even for zero 
vertical emittance we never have full coherence in the vertical direction. This 
unexpected result is due to the influence of the horizontal emittance on the 
vertical coherence properties of the photon beam. Because of this, the degree 
of coherence changes between zero and unity within the diffraction angle. We 
also studied the near field zone. When one is interested in the evolution of the 
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degree of coherence along the beamline back up to the exit of the undulator 
the situation becomes much more complicated with respect to the far zone 
case, as the observation distance is one of the problem parameters. There are 
many more asymptotic situations which can be studied, and a large part of 
our paper is devoted to the calculation of these asymptotic situations. We 
provided approximate estimations for the vertical coherence length that are 
valid from the far zone and back, up to the exit of the undulator in the case 
when either the vertical Fresnel number or the vertical electron beam diver- 
gence are much larger than unity. These can be used at the stage of planning 
experiments. 

It should be noted that, throughout this work, we did not discuss the accuracy 
of the approximation of small and large parameters. In order to do so, one 
needs to develop a perturbation theory for each asymptotic case studied here, 
which would considerably increase the size of this paper. As a result we leave 
this issue for future work. 

Finally, we selected an application to show the power of our approach. We 
discussed how the transverse coherence properties of an X-ray beam can be 
manipulated to obtain a convenient coherent spot-size on the sample with the 
help of a simple vertical slit; this invention was predicted almost entirely on 
the basis of theoretical ideas of rather complex and abstract nature discussed 
in the previous parts of the paper. 
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Appendix A: Random phasor sum 



The field of thermal light can be regarded as a sum of a great many indepen- 
dent contributions. The complex envelope of polarized thermal light at fixed 
time and a fixed point in space is a sum of a very large number of complex 
phasors 

N 

£(f,t) = ]Ta fc e* , (157) 

k=l 



where N is the number of radiating atoms. Statistical properties of elementary 
phasors that are generally satisfied in thermal light problems of interest are 
as follows: 

a) The amplitudes a& = Re(a>k) and the phases ipk are statistically indepen- 
dent of each other and of the amplitudes and phases of all other elementary 
phases for different values of k. 

b) The random variables are identically distributed for all k with mean 
value < a > and second moment < a 2 >. 

c) The phases ipk are uniformly distributed over the interval (0, 2n). 

The reader can find in H that when assumptions from a) to c) are satisfied, 
the real (Re(E)) and imaginary (Im(E)) parts of the field are distributed in 
accordance with the Gaussian law in the limit for iV — > oo, so that 



[Re(E)f 



(ME)) 
[Im(E)f 



(Im(E)) 



a 



p(Re{E),Im(E)) 



2tt(j'- 



exp 



{Re(E)Im(E)) = 
[Re(E)] 2 + [Im{E)] 2 



2a 2 



;i58) 



where p(Re(E), Im(E)) is the joint probability density function. 

In Section 2.1 we discussed statistical properties of Synchrotron Radiation and 
we were led to assumptions 1), 2) and 3) which are weaker than a), b) and c). 
Here we will demonstrate that assumptions from a) to c) can be relaxed to 
assumptions from 1) to 3) without changes in results. We will derive results 
valid when the amplitudes ol\~ are complex = \ak \ exp (i(j)k)- 

After denoting with r and i the real and imaginary parts of the fields and 
after substituting notation (Q) with Q we first demonstrate that f = i = 0. 
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We have straightforwardly 



1 N ( \ 

r = —^2 \\ a k \ cos <f) k cos ip k - \ a k \ sm<f) k smip k 1=0, 

1 N / \ 

i = — (J a k I cos 0fe sin^fe + | a. k | sin0 fc cos^fcj = , (159) 

because all averages over trigonometric functions are zero. 

Second, we demonstrate that r 2 = i 2 = a 2 /2. Again, direct calculation shows 



N 



N 



k,n=l 





On | 


COS <p k COS (f) n COS Tpk COS l]) n 


+ ttfc 




sin (p k sin <p n sin ip k sin ip n 


- OL k 


\ O n 


cos (p k sin (fr n cos ip k sin ip n 



-\a k \\a n \ sin0 fe cos0„ sin tp k sin ijj n J 

— £ I «fc I 2 (cos 2 ^ fc + sin 2 ^fc) = y ■ (160) 
fc=i 



Moreover it is easy to see that i 2 = r 2 . 
Finally, we show that ri — 0. In fact 



_ 1 

ri — — 



a k 




COS <j) k COS <j) n 


cos i\) k sin i\) n 


ot k | 




cos 4> k sin <f) n 


COS 1p k COS 1p n 




1 a n 


sin 0^ cos 4> n 


sin ip k sin 



2AT 



a k \ \ a n \ sin0 fc sin0 n sin^ fc cos^ n 



£ |a fc | 2 (cos0 fe sin0 fc - sin0 fc cos0 fc ) = . (161) 



As a result we have that real and imaginary parts have zero means, equal 
variances and are uncorrelated. Use of the central limit theorem allows to 
conclude that the resulting phasor sum is a circular complex Gaussian random 
variable. 
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Appendix B: A useful transformation of the expression for the un- 
dulator radiation field 



We start reporting here, for convenience, Eq. (24), that represents the field 
(in normalized units) produced by a particle with offset and deflection at any 
distance z ^ 1/2 from the exit of the undulator, where the undulator center 
is taken at z Q = 0: 



1/2 

= z a J dz 

-1/2 



Z n - Z 



- exp < i 



r± -l- fjz' 



21 



{162) 



In this Appendix we show that E s ± as reported in Eq. (162) may be described 

as 



1/2 



f z D dz' 
= J i _ -r, exp \ 



Z n - Z' 



-1/2 

where $[/ is given by 



$c/ + Cz' + 



Z n Z 



2(z Q - z<) 





y 


2- 


> 


M 


- V 




> 




J 




> 



> (163) 



u 



Qy ~ 



Zo 

2 



(164) 



This means that E s ± is of the form 



E s± C, z , 6 - (l/z ) -r})= exp(i^u)S C, z Q , - (l/z ) - fj) . (165) 



Let us introduce in this Appendix, for simplicity of notation, £ = fj_ — I and 

— * 

4> = i/z Q . It is easy to rewrite the phase in the integrand of Eq. (162), which 
we denote with $ T as 



^2 



<$>T = Cz'+\z> + \ 



z - 2(f) ■ fjz' + 



■i] z 

Zo 



12 



oo f% 



1 + 



(166) 



n=l 



Further algebraic manipulation of Eq. (166) yields 
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->2 



,2 oo 



+0 4£ f f -(167) 



n=2 



n=l \ Z °, 



It is easy to recognize $j/ in the second term on the right hand side of Eq. 
(167). Furthermore, the fourth term on the right hand side of Eq. (167) can 
be further manipulated, leading to 



$ T = Cz! + $u + — ( (f) - fj 



*\ ^ ( z' 



(168) 



n=2 



that is 



$ T = C7i' + $[/ + — 



2 oo / 2'\ US 



n=2 \ Z ° 



(169) 



or 



2 Z 



2(4 - 5') 



(170) 



Therefore, since 



— fj = 9 17 



(171) 



we have 



1/2 

/■ z dz' 

J Z n — Z' 



exp < 



-1/2 



ZnZ' 



2(4 - 5') 



^ i 

9---fj 

Z n 



(172) 



that is Eq. (163), quantum erat demonstrandum. 
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Appendix C: Autocorrelation function for undulator sources 



In this Appendix we demonstrate the validity of Eq. (114). Having defined 



oo oo 



f(z Q , A9 y ) = J 4 V f dkS* [z , <g + {4> y - A9 y ) 2 ]s[z , <g + (0, + Aft 



-oo — oo 



and 



oo oo 



2tt 2 



d(f)y / d(f) x sinc 



-oo — oo 



4>l + (4>y ~ AOyf 



sine 



we have want to demonstrate that Eq. (114) holds, that is 
f(z ,A9 y )=P(A9 y ) . 



(175) 



The proof is based on the autocorrelation theorem, which states that if the 
(two-dimensional) Fourier Transform of a function w(x,y) with respect to 
variables a x and a y is indicated by w(a x ,a y ), then the Fourier transform of 
the two-dimensional autocorrelation function of w(x,y) with respect to the 
same variables a x and a y is given by \w(a x , a y )\ 2 . In formulas, after definition 
of the autocorrelation function 



oo oo 



A[w](x,y)= J drj J d^w(rj + x,£ + y)w*(rj,£) , 



(176) 



— oo — oo 



which is equivalent to 

oo 

A[w](x,y)= [ d V I d{w{ri + x/2,S + y/2)w*{ri-x/2,S-y/2) , (177) 



oo oo 



— oo — oo 



the autocorrelation theorem states that 



oo oo 



J dx J dyexp [i(a x x + a y y)]A[w](x , y) = \w(a x ,a y )\ 2 . (178) 



-oo — oo 



First we extend the definition of / 
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oo oo 



f(z , A9' x , A9' y ) = JL j d&y j d4> x S* [zo, (4 - A^/2) 2 + (0, - A6' y /2) 2 



-oo — oc 



S o ,(^ + A0;/2) 2 + (0, + A0'/2)' 



(179) 



where we changed variables from A9 XjV to = 2A0 x>y Then we can apply 

the autocorrelation theorem in Eq. (178) to the function / thus obtaining the 
following relation: 



oo oo 

dA9' x J dA^exp [ l (a^ + a^;)]/(i,AtA^ 

X. — oo 

oo oo 

12 i ?2 



27T 2 



J d(f) x J d<j> y exp [i(a x (j> x + ai^)]^ z , 2 + 2 



Hoo — oo 



,(180) 



where a XjV are now conjugated variables with respect to the angles <f) x>y on 
which S depends. We will denote with S the two-dimensional Fourier Trans- 
form of S, that is: 



oo oo 



S(a x ,a y )= J d<j) x j c%exp [i(a x (f) x + a y <j) y )}S z a ,(j) 2 x + <j) 2 y . (181) 



-oo — oo 



The relation between the function S and the undulator field is given by Eq. 
(31), and one has 



oo oo 



S(a x , a y ) = J d(j) x J d(j) y exp [i(a x (j> x + a y (f) y )] exp [-i(4> 2 x + 4> 2 y )z /2\ 



-oo — oo 



XE 



Zo,</>l + <P 2 y 



;i82) 



After definition of a XjV = a X)V /z one can write S, as a function of a X:V instead 
of a Xj y. Then, one can switch to the new integration variables x = 4> x z and 
y = (f) y z to obtain: 



S(a x ,a y ) = — 



dx j dy exp [i(a x x + a y y)} exp [— i(x + y)/(2z )\ 

xE s± z Q ,x 2 + y 2 



(183) 



where the expression for E s j_[z , x 2 + y 2 } is given in Eq. (24). Now we have to 
calculate the Fourier transform of the product of two factors: exp [— i(x 2 + y 2 ) / (2z )\ 
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and E ±s . Let us look for the Fourier transform of each factor. 
A direct calculation shows that 



oo oo 



J dx J dyexp [i(a x x + a y y)]E s ± (z , x 2 + y 2 ) 



-oo — oo 



2m z o exp 



— i 



[a 2 x + a 2 y )z 



sine 



a 2 x + a 2 y 



(184) 



Second, let us deal with the Fourier transform of exp [— i(x 2 + y 2 )/(2z )\. 



oo oo 



J dx j dyexp [i(a x x + a y y)} exp [— i(x 2 + y 2 )/(2z ) 



-oo — oc 



= -4iz Q exp [i(a 2 x + a 2 y )z Q /2] . (185) 

Since the Fourier transform of a product is equal to the convolution of the 
Fourier transforms of the factors we have 



oo oo 



S(a x , a y ) = 8ir J du J dw exp {i[(a x — u) 2 + (a y — w) 2 ]z Q /2} 

— oo — oo 

x exp \-i(u 2 + w 2 )z /2]smc[(u 2 + w 2 )/A] . (186) 

and therefore we have 



\S(a x ,a y )\ = 8tt 



oo oo 



du / dwexp [i(a x u + a y w)z ]sinc[(u + w )/4][ L87) 



-oo — oo 



Going back to old variables a XjV we obtain 



\S(a x ,a y )\ = 87r 



oo oo 



J du J dwexp [i(a x u + a y w)]smc{(u 2 + w 2 )/A] 



-co — oo 



.(188) 



which is independent of z . As a result Eq. (180) is also independent on z Q , 
i.e. the Fourier transform of / is independent of z Q . 

Now, on the one hand in the limit for z Q — ► oo the function / transforms 
into [3, because the S functions in / tend asymptotically to the sine functions 
in f3. On the other hand, if the Fourier transform of / is independent of z a , 
also / is independent of z Q . As a result it can only be f(A9' y ) = f3(A9' y ), and 
f(A0y) = /3(A0y) that is Eq. (175) holds, quantum erat demonstrandum. 
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Note that, based on the autocorrelation theorem, it is also possible to give 
an analytic expression for the Fourier transform of (3. After definition of f3 — 
(3(z Q , A9' x , A9' y ) as in Eq. (179), application of the autocorrelation theorem 
simply states that the two-dimensional Fourier Transform of (3, that will be 
indicated with f3 can be written as 



(3(a x , a y 



1 



2tt 2 



d<p x / d(j) y exp [i(a x (f) x + a y (f) y )}smc 



+ 



.(189) 



Introducing a 2 = a 2 x + a 2 and representing the two-dimensional Fourier Trans- 
form of (3 in terms of Fourier-Bessel transform we obtain 



/3(a) 



2tt 2 



2tt J dr rJ (ra)sinc (^^j 



= 2 



7T - 2Si(a 2 ) 



(190) 



where Si indicates the sine integral function. 

Finally, one can get back a simpler representation of the function (3 in terms of 
a one-dimensional integration simply performing an anti Fourier-Bessel trans- 
form: 

1 °° 

(3{A9') = — { da aJ (aA9')(3(a) , (191) 
2n J 



where A# 



A9J + A9 y 2 . For A9 X = we obtain 



if ( AO \ 

(3(A0y) = -J daaJoi a^- J [tt - 2Si(a 



(192) 
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